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Propagation  of  Intense  Short  Laser  Pulses 
in  the  Atmosphere 


P.  Sprangle,  J.  R.  Penano  and  B.  Hafizi* 
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Naval  Research  Laboratory 
Washington  D.C.  20375 

^Present  address:  Icarus  Research,  Inc.,  P.O.  Box  30780,  Bethesda,  MD  20824-0780. 

Abstract 

The  propagation  of  short,  intense  laser  pulses  in  the  atmosphere  may  have  a  number  of 
important  applications  in  the  areas  of  active  and  passive  remote  sensing,  electronic  countermeasures 
and  induced  electric  discharges  (artificial  lightning).  For  example,  localized  ultraviolet  radiation 
generated  at  a  remote  distance  by  laser  ionization  can  provide  a  source  for  active  fluorescence 
spectroscopy  of  biological  and  chemical  agents  in  the  atmosphere.  Pulses  of  intense,  directed  white 
light  may  also  find  applications  in  the  areas  of  hyperspectral  imaging  and  differential  absorption 
spectroscopy.  In  this  paper  the  propagation  of  short,  intense  laser  pulses  in  the  atmosphere  is 
investigated  theoretically  and  numerically.  A  set  of  3D,  nonlinear  propagation  equations  is  derived, 
analyzed,  and  solved  numerically.  The  propagation  equations  include  the  effects  of  dispersion, 
nonlinear  self-focusing,  stimulated  molecular  Raman  scattering,  multi-photon  and  tunneling 
ionization,  energy  depletion  due  to  ionization,  relativistic  focusing,  and  ponderomotively  excited 
plasma  wakefields.  The  instantaneous  frequency  spread  along  a  laser  pulse  in  air,  which  develops 
due  to  various  nonlinear  effects,  is  analyzed  and  discussed.  Coupled  equations  for  the  power,  spot 
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size  and  electron  density  are  derived  for  an  intense  ionizing  laser  pulse.  An  approximate 
equilibrium  of  a  single  optical  filament,  which  involves  a  balancing  between  diffraction,  nonlinear 
self-focusing,  and  plasma  defocusing,  is  obtained  and  shown  to  require  a  specific  distribution  of 
power  along  the  filament.  These  equations  indicate  that  in  the  presence  of  ionization  a  self-guided 
optical  filament  is  not  realizable.  A  method  for  generating  a  remote  spark  in  the  atmosphere  is 
proposed.  This  method  utilizes  the  dispersive  and  nonlinear  properties  of  air  to  cause  a  low 
intensity  chirped  laser  pulse  to  compress  both  longitudinally  and  transversely.  For  optimally  chosen 
parameters,  the  transverse  and  longitudinal  focal  lengths  can  be  made  to  spatially  coincide,  resulting 
in  rapid  intensity  increase,  ionization,  and  white  light  generation  in  a  localized  region  far  from  the 
source.  Coupled  equations  for  the  laser  spot  size  and  pulse  duration  are  derived  that  can  describe 
the  focusing  and  compression  process  in  the  low  intensity  regime.  More  general  examples 
involving  beam  focusing,  compression,  ionization,  and  white  light  generation  near  the  focal  region 
are  studied  by  numerically  solving  the  full  set  of  3D,  nonlinear  propagation  equations. 


VI 


I.  INTRODUCTION 


Experiments  using  ultra  short  (-100  fsec),  high  intensity  (>  10^^  W/cm^)  laser  pulses 
appear  to  have  demonstrated  long-distance  atmospheric  propagation,  air  breakdown, 
filamentation,  and  white  light  generation  [1-14],  Intense,  directed  white  light  pulses  have  been 
generated  and  backscattered  from  atmospheric  aerosols  up  to  altitudes  of -15  km  [4],  The 
generation  of  pulsed  THz  radiation  in  plasma  channels  formed  by  femtosecond  pulses  has  also 
been  observed  and  analyzed  [15,16],  Although  many  of  the  observations  can  not  be  completely 
explained,  the  experimental,  theoretical,  and  numerical  results  obtained  to  date  indicate  potential 
applications  for  both  passive  and  active  remote  sensing  [4,  17,  18],  and  induced  electric 
discharges  [10-14],  among  others.  To  achieve  these  potential  applications  it  is  necessary  to  have 
a  comprehensive  and  quantitative  understanding  of  the  physical  mechanisms  that  govern  the 
propagation  of  intense,  short  laser  pulses  in  air. 

The  propagation  of  intense,  short  laser  pulses  in  the  atmosphere  involves  a  variety  of 
diverse  linear  and  nonlinear  optical  processes.  The  combined  effects  of  diffraction,  nonlinear 
self-focusing,  ionization,  and  plasma  defocusing  play  an  important  role  in  the  propagation  of 
laser  and  plasma  filaments  [1,  5,  6,  19-26].  In  addition,  nonlinear  bound  electron  effects, 
stimulated  Raman  scattering,  and  plasma  formation  contribute  to  considerable  spectral 
broadening  and  white  light  generation  by  the  laser  pulse  [8,  9,  23,  27-32]. 

The  physics  governing  the  atmospheric  propagation  of  short  intense  laser  pulses  can  be 
very  different  from  that  of  long  laser  pulses.  For  example,  the  Raman  instability  associated  with 
the  excitation  of  molecular  rotational  modes,  which  can  disrupt  the  long  distance  propagation  of 
long  (>  nsec)  pulses  [28],  may  not  be  as  disruptive  for  laser  pulses  that  are  shorter  than  the 
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characteristic  period  of  the  rotational  mode  (~  psec).  In  addition,  experiments  and  theory 
indicate  that  the  nonlinear  refractive  index  of  air  is  a  fiinction  of  the  laser  pulse  length;  e.g.,  for  a 
-100  fsec  pulse,  it  is  observed  that  the  effective  nonlinear  refractive  index  can  be  several  times 
smaller  than  for  a  longer  (>  psec)  pulse  [1],  Also,  because  of  their  large  spectral  content,  short 
laser  pulses  are  more  affected  by  dispersion.  Finally,  the  atmospheric  propagation  of  intense, 
short  laser  pulse  trains  generated  by,  for  example,  an  RF  linac  driven  free  electron  laser  [33], 
may  result  in  sufficient  spectral  broadening  to  affect  the  laser  absorption  rate.  That  is,  the 
broadened  laser  pulse  spectrum,  rather  than  lying  between  individual  absorption  lines,  may 
overlap  some  of  the  lines.  This  could  affect  the  thermal  blooming  process  which  is  a  sensitive 
function  of  the  absorption  rate. 

In  this  paper  we  derive,  analyze  and  numerically  solve  a  system  of  3D,  nonlinear 
equations  for  atmospheric  laser  pulse  propagation.  The  model  includes  diffraction,  group 
velocity  and  higher  order  dispersion,  stimulated  molecular  Raman  scattering,  photoionization, 
nonlinear  bound  electron  effects,  ionization  energy  depletion,  and  propagation  in  a  spatially 
varying  atmosphere.  The  propagation  equations  are  used  to  analyze  a  number  of  physical 
processes,  such  as  optical/plasma  filamentation,  pulse  compression,  nonlinear  focusing  and  white 
light  generation,  see  Fig.  1.  A  coupled  set  of  equations  for  the  laser  power,  spot  size,  and  electron 
density  is  derived.  At  sufficiently  low  intensities,  where  ionizaton  is  neglegible,  a  neccessary  condition 
for  equilibrium  of  a  laser  filament  is  derived.  It  is  shown  that  at  high  intensities  laser  power  depletion  due 
to  ionization  implies  the  absence  of  a  matched  beam  solution,  although  extended  propagation  is  possible 
provided  the  ionization  rate  is  sufficiently  small.  For  laser  intensities  sufficiently  low  that  ionization 
effects  and  stimulated  Raman  scattering  can  be  neglected,  a  set  of  coupled  equations  for  the  laser 
spot  size  and  pulse  duration  are  also  derived.  The  coupled  equations  describe  nonlinear  self- 
focusing,  compression  and  spreading  of  chirped  pulses  in  a  spatially  varying  atmosphere. 
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Simulations  based  on  the  3D  numerical  solution  of  the  general  propagation  equations  are  used  to 
study  highly  nonlinear  propagation  in  the  presence  of  plasma  generation  and  Raman  scattering. 

A  process  by  which  a  laser  pulse  can  remotely  ionize  a  localized  region  of  the 
atmosphere  is  studied  for  possible  remote  sensing  applications.  By  introducing  a  negative 
frequency  chirp  on  a  relatively  long  laser  pulse  (>  psec)  the  pulse  can  undergo  longitudinal 
compression,  due  to  linear  group  velocity  dispersion.  In  addition,  transverse  self-focusing  of  the 
pulse  takes  place  due  to  atmospheric  nonlinearities.  For  a  properly  chosen  set  of  parameters,  the 
focal  distances  for  longitudinal  compression  and  transverse  focusing  can  be  made  to  coincide, 
resulting  in  a  significant  intensity  increase  over  a  relatively  localized  region.  The  compressed 
and  focused  laser  pulse  can  ionize  a  local  region  (~  Im  in  extent)  of  the  atmosphere  kilometers 
away  from  the  source.  The  localized  spark  can  generate  ultra-violet  radiation  through 
recombination.  Since  many  biological  and  chemical  agents  will  fluoresce  in  the  optical  regime 
when  illuminated  with  ultra  violet  radiation,  the  recombination  radiation  can  be  used  as  a  source 
for  atmospheric  fluorescence  spectroscopy.  For  this  application,  the  ultra-violet  radiation  must 
be  generated  locally  since  it  is  highly  absorbed  in  the  atmosphere. 

This  paper  is  organized  as  follows.  In  Sec.  II  the  general  nonlinear  3D  propagation 
equations  are  derived.  In  Sec.  Ill  photoionization  processes,  filamentation  in  neutral  air  and 
white  light  generation  (including  the  effects  of  photoionization  and  Raman  scattering)  are 
discussed.  Using  the  source-dependent  expansion  method,  coupled  equations  for  the  laser 
power,  spot  size  and  electron  density  are  derived,  and  an  approximate  equilibrium  obtained  in 
Sec.  IV.  Section  V  presents  a  discussion  of  compression  and  focusing  of  a  laser  pulse  in  the 
atmosphere,  leading  to  the  generation  of  a  spark  at  a  remote  location.  A  set  of  coupled  equations 
for  the  laser  spot  size  and  pulse  duration  are  derived  from  the  propagation  equations  in  the  low 
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intensity  regime.  Using  the  simplified  coupled  equations,  conditions  are  derived  for  optimal 
compression  and  focusing  in  the  low  intensity  limit.  Simulations  based  on  numerical  solution  of 
the  general  propagation  equations  on  a  3D  Cartesian  grid  are  used  to  model  high  intensity 
atmospheric  propagation.  Laser  pulse  propagation  in  a  spatially  varying  atmosphere  is  also 
considered.  Finally  a  summary  is  presented  in  Sec.  VT.  Appendix  A  contains  the  derivation  of 
the  linear  source  terms  of  the  propagation  equation  while  appendix  B  contains  the  derivation  of 
the  nonlinear  source  terms. 


IL  GENERAL  NONLINEAR  PROPAGATION  EQUATION 

In  this  section  we  derive  a  general  nonlinear  3D  equation  describing  the  propagation  of 
an  intense  laser  pulse  in  air.  The  equation  incorporates  the  effects  of  diffraction,  dispersion, 
ionization,  pulse  energy  depletion  due  to  ionization,  stimulated  molecular  Raman  scattering, 
nonlinearities  associated  with  bound  electrons,  spatial  inhomogeneity  in  air  density,  plasma 
Wakefields  and  relativistic  electron  motion. 


i)  Wave  Equation 

The  starting  point  is  the  wave  equation  for  the  laser  electric  field  E(/',  /),  given  by 


where  Vj  is  the  transverse  Laplacian  operator  and  z  is  the  coordinate  in  the  direction  of 
propagation.  The  quantities  S/,  and  S denote  source  terms  which  are  respectively  linear  and 
nonlinear  in  the  laser  electric  field. 


4 


The  laser  electric  field,  y,  z,  t) ,  linear  source  term,  S/,  (x,  y,  z,  t) ,  and  nonlinear 
source  term,  S]yi{x,y,z,t) ,  are  written  in  terms  of  complex  amplitudes,  A{x,y,z,t), 

Sl  (x,  y,  z,  t)  and  S (x,  y,  z,  t)  and  a  rapidly  varying  phase,  y  (z,  t) ,  that  is 

E(x,  y,  z,  0  =A(x,  y,  z,  t)  exp(z  y/{z,  11  +  c.c. ,  (2a) 

Si{x,y,zJ)^Si{x,y,z,t)Qxy{iy/{zA))txll  +  c.c.,  (2b) 

S NL  0  =‘S'iw:  0 exp(/ t))  tx/2  +  c.c. ,  (2c) 

where  y/{z,  i)  =  kg  z  -  cOot  is  the  phase,  kg  is  the  carrier  wavenumber,  (Og  is  the  carrier 
frequency,  e^.  is  a  transverse  unit  vector  in  the  direction  of  polarization,  and  c.c.  denotes  the 

complex  conjugate.  Substituting  the  field  and  source  representations  given  by  Eqs.  (2)  into  Eq. 
(1)  yields. 


~  ^ - f"  +  2i  kg 

c^ 


^AOg 
+  2/  — 


c  dt 


^(r,0  =  +  5yii{Y,{), 


(3) 

where  the  rapidly  varying  phase  factor  has  been  cancelled  from  both  sides  of  the  equation. 
Although  the  atmospheric  density  is  spatially  varying,  the  wavenumber  is  taken  to  be  constant 
since  the  maximum  change  in  the  linear  refractive  index,  i.e.,  fractional  change  in  wavenumber, 

from  sea  level  to  vacuum  is  <  10“'^ .  The  expressions  for  the  linear  and  nonlinear  source 
amplitudes  are  derived  in  appendices  A  and  B,  respectively, 
ii)  Linear  Source  Terms 

The  linear  source  amplitude  can  be  expressed  as 


Yj'  cc,ir)  6);' 


^=0 


d'AjrJ) 

dt' 


(4) 
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where  I  -  0,\,2,  ■■■ .  The  unitless  dispersion  coetFicients  (r)  in  Eq.  (4)  are  given  by 


ai  =  - 


(O 


i-2 

0 


(5) 


where  P{co)  ^  {co I c){  \  +  Ak XLi.O))\' ^  =  (<y/c)«o(<w),  XLi(o)  is  the  linear  susceptibility  of 
bound  electrons  and  (co)  is  the  refractive  index. 

Hi)  Nonlinear  Source  Terms 

The  nonlinear  source  amplitude  is  due  to  a  number  of  effects  and  can  be  written  as 

~~  ^ bound  ^ Raman  ^plasma  ^wake  ^rel  ^wn->  (^) 

where  the  individual  contributions  are  described  as  follows.  The  nonlinear  contribution  from 
bound  electrons,  Kerr  effect,  is  given  by 

=  (7) 

where  is  the  electronic,  contribution  to  the  nonlinear  refractive  index.  The  nonlinear  index 

defines  a  nonlinear  self-focusing  power  [34-36],  Ppn  = 

The  source  term  due  to  stimulated  molecular  Raman  scattering  is  given  by 

2 

^Ramani^>t)^  -  4;r 0(O^(r,O,  (8) 

c 

where  xi  the  linear  susceptibility  and  the  unitless  Raman  oscillator  function  Q(i)  is 
determined  by  solving  Eqs.  (B9,  BIO).  The  Raman  source  term  can  also  contribute  to  the  third 
order  polarization  field. 

The  plasma  source  term  is  given  by 
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^ plas)7tai^^0  ~  2 


2  ^  \ 

1-/-^ 
COq 


Air,t), 


(9) 


(  2 

where  o)p (r, 0  =  \4;t  q  (r, /) /m j  is  the  plasma  frequency,  is  the  plasma  density 


generated  by  ionization,  and  is  the  electron-neutral  collision  frequency.  Ionization  results  in 

a  plasma  column  which  is  localized  to  the  laser  axis.  The  plasma  column  causes  a  local  decrease 
in  the  refractive  index  which  can  defocus  the  laser  pulse.  The  term  proportional  to  the  electron 
collision  frequency  is  responsible  for  the  collisional  absorption  of  laser  energy,  i.e.,  inverse 
bremsstrahlung. 

The  source  term  Sy^ake  is  due  to  the  possible  generation  of  plasma  waves  and  is  given  by 


^wake  (*”>  0  ~ 


(Op{r,t)  dn, 


A{r,t), 


(10) 


where  Sn^  represents  a  plasma  density  perturbation  driven  by  the  ponderomotive  force  of  the 

laser  pulse,  i.e.,  a  plasma  wakefield  [37].  The  density  perturbation,  determined  by  Eq.  (B19) 
together  with  Gauss's  equation,  results  in  a  modulation  of  the  plasma  density  at  the  plasma 
frequency. 

The  term  is  due  to  relativistic  effects  arising  from  quiver  motion  of  plasma  electrons 
in  the  field  of  the  laser  and  is  given  by 


Srelirj) 


4c^ 


q\A{r,t)\ 

mca>Q 


(11) 


This  relativistic  source  term  defines  a  critical  self-focusing  power  due  to  plasma  [36], 

Pplasma  =  2c{ql r^Y  «o{<»o  IcOp'Y  where  I mc^  is  the  classical  electron  radius.  The 


total  nonlinear  self  focusing  power  consists  of  contributions  from  both  Pjvi  and  P plasma  and  is 
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given  by  P plasma  KPnl  +  Pplasma  )  as  shown  in  Ref.  36.  Typically  Ppiasma  »  Pnl^  so 

that  the  nonlinear  self  focusing  power  is  due  to  bound  electrons  and  is  equal  to  Ppn  ■ 

Finally,  the  term  describing  the  depletion  of  laser  energy  due  to  ionization  is  given  by 

=  (12) 

where  U ion  is  the  characteristic  ionization  energy.  For  example,  the  ionization  energy  for  O2  is 
12.1  eV  while  for  N2  it  is  15.6  eV. 

iv)  Full  Nonlinear  Three-Dimensional  Propagation  Equation 

Substituting  Eqs.  (4)-(12)  into  Eq.  (3)  results  in  the  following  nonlinear  propagation 
equation  for  the  laser  envelope. 


\-i- 


co 


O  J 


1.1  k - 1 - — 

dz  dz^ 


+  2/  (1  -«!  /  2)^  -  (1  -^2 )-rT2 

c  dct  dc  t 

1  . 1 2  ,  \  a\2  _ 

",  ■  4;rc  I  '  ' 


+  ioc2 


(Oo  dc^P 


]A{r,t)  = 


act  c. 


(13) 

2  2  2  2 

where  the  summation  in  Eq.  (4)  has  been  limited  to  ^<3  and  AK  =  (l-or^  )<»(,  /c  —  ko  . 

It  proves  useful  to  transform  the  independent  variables  from  z,  t  to  z,  r,  where 
T  =  t-  zNg  and  v^  will  be  set  equal  to  the  linear  group  velocity  of  the  pulse.  In  terms  of  the 


new  variables  the  derivatives  transform  as,  d! dt  — >  5/ dr  and  d! dz  d! dz  -  v g  d! dr . 

Under  this  transformation,  Eq.  (13)  becomes 
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[  - 


® , 


+  2i 


2  2 


Anc 


-Mf 


2  _  P 

1-/- 

c"  1 

AD  a 

C  dCT  ^ 

+  2ik, 


a  2  5'  a' 


+ 


dz  dzdcx  dz‘ 


a^  .  c  d- 


a  2  2 

dc  T 


+  la^ 


OJo  dc^T^ 


]A{x,y,z,T) 


4m^  cOq 


,2 


A 


Uio„  ^  +  ^^^XLQit) 

rig  OCT 


A(x,y,z,T), 

(14) 


where  A/2  -  (l-aj/ 2)(0q  -  Pg^  ckg  and  Pg  =  Wg  Ic  . 

The  wavenumber  k^  and  group  velocity  Vg  appearing  in  Eq.  (14)  are  as  yet  unspecified. 

It  is  convenient  to  choose  them  so  that  the  form  of  the  propagation  equation  is  simplified. 
Choosing  AK  =  0  and  A/2  =  0  defines  the  carrier  wavenumber  and  linear  group  velocity 

respectively  as  A:,,  =  (l-a,,)  6)^,  / c  =  n^cOg  / c ,  and 

Vg  =  cn^  /(I- ai)  =c /(«£,  +  cOg, dn^  I d co q)  .  T aking  Ml  =  A/2  =  0 ,  the  propagation  equation 
simplifies  to 


[  V' 

L  „2 

^  V 


1  -i 


V, 


CO 


a  2  a"  a" 

+  2/A:, - + 


oj 


°  dz  p  dzdcx  dz^ 


KPl  o  +  ^^3  — - ^  ]A(x,y,z,T) 


dc^T^ 


2  2 
4rrc 


Mf 


CO, 


+- 


4m^  a)^ 


MP  - 


\ 


dcT  r.^ 


J 


A{x,y,z,T) 

(15) 


9 


Equation  (15)  describes  the  3D  evolution  of  the  complex  laser  field  amplitude,  A(x,  y,  z,  r) .  The 
self-consistent  model  employed  here  involves  the  solution  ofEq.  (15)  along  with  equations  that 
describe  the  response  of  the  medium  (air)  to  the  laser  field.  In  Eq.(15),  the  linear  dispersion 
response  is  obtained  from  Eq.  (A  16),  while  the  nonlinear  bound  electron  response  is  given  by  Eq. 
(B6).  The  stimulated  Raman  response  is  obtained  from  Eqs.  (B9),  (BIO)  and  (Bll).  Plasma 
effects,  Wakefields  and  relativistic  effects  are  given  in  Eq.  (B20)  while  Eq.  (B26)  provides  the 
expression  for  pulse  energy  depletion  due  to  ionization.  Finally,  photoionization  rates  are  given 
in  the  following  section. 

A  3D  numerical  simulation  based  on  solving  Eqs.  (15)  together  with  the  medium  reponse 
has  been  developed  which  places  the  laser  pulse  on  a  Cartesian  (x,y,x)  grid,  allowing  for  the 
modeling  of  asymmetric  pulse  shapes  and  laser  filamentation.  The  laser  pulse  is  advanced  in  z 
according  to  Eq.  (15)  using  a  split-step  method  [38]  in  which  the  linear  terms  are  advanced  in 
Fourier  space,  while  the  nonlinear  terms  are  handled  in  coordinate  space.  The  equations 
describing  ionization,  wakefield  generation  and  Raman  scattering  are  solved  at  each  z-step  by  a 
4*  order  Runga-Kutta  integration.  To  facilitate  computation  of  the  plasma  wakefield,  the  term 

c^VxVx  has  been  neglected  and  the  approximation  V  «  V_l  -  (e^  Iv  g){dl  dr)  has  been 
made  in  Eq.  (B19)  and  in  using  Gauss’s  equation. 

III.  IONIZATION,  FILAMENTATION  AND  WHITE  LIGHT  GENERATION 

In  the  following  subsections  we  use  the  theoretical  model  presented  in  the  previous 
section  to  analyze  ionization,  filamentation,  and  spectral  broadening  of  short,  intense  laser  pulses 
in  air.  The  analyses  of  optical  filament  propagation  and  white  light  generation  consist  of 
substituting  a  self-similar  form  for  the  solution  of  the  complex  amplitude  into  a  reduced  version 
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of  Eq.  (15)  and  obtaining  equations  for  such  quantities  as  the  laser  spot  size,  phase,  curvature, 
and  instantaneous  frequency.  While  these  analyses  are  not  meant  to  be  rigorous,  they  can 
provide  some  quantitative  understanding  of  the  processes  considered. 

i)  Photo-Ionization 

The  free  electron  density  in  air  can  change  because  of  ionization,  recombination  and 
attachment  processes.  The  rate  equation  for  electron  density  «e  is 

^  =  Wn„  -Tin.-Pr-nl,  (16) 

ot 

where  is  the  neutral  gas  density,  W  is  the  photoionization  rate,  77  is  the  electron-attachment 
rate  coefficient  and  is  the  recombination  coefficient.  Empirical  relationships  for  the 
attachment  and  recombination  rates  in  air  are  available  [39].  The  recombination  time  is 

Tf.  =l/(y!?r/7g)  ,  where  ~  3.7x10”^  cm^/sec  [14],  As  an  example,  for  a  plasma  with 

density  of  =  lO^^cm”^ ,  the  recombination  time  is  Ty  «  2.7 nsec .  The  attachment  time  is 
expressible  as  =  1  /  7 ,  and  for  typical  atmospheric  parameters  and  laser  intensities  considered 
here  r  4  is  on  the  order  of  1  psec.  Since  the  laser  pulses  of  interest  here  have  a  duration  of  1 
psec  or  less,  recombination  and  attachment  processes  play  a  negligible  role  in  the  propagation  of 
a  single  pulse. 

For  short  laser  pulses,  free  electrons  are  generated  by  multi-photon  and  tunneling 
processes;  avalanche  ionization  is  not  significant.  The  principal  constituents  of  air  are  nitrogen 
and  oxygen  and  hence  the  total  photoionization  rate  can  be  written  approximately  as 
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W  =  +  0.2  ,  corresponding  to  the  proportion  of  A^2  and  O2  molecules  in  normal 

atmosphere.  For  either  one  of  these  species  the  ionization  rate  takes  different  forms  according  to 
the  value  of  the  Keldysh  parameter  xk  t^O],  i.e., 


rK  =  2.3lxW 


;i^[pm]/[W/cm^ 


where  UiQ„  is  the  ionization  energy.  The  value  of  the  Keldysh  parameter  identifies  the  multi¬ 
photon  (f ^  » 1 )  and  the  tunneling  {yk  «  0  regimes.  In  the  multi-photon  regime  the 
ionization  rate  is  [41] 


2n(Oo  I{r,z,T) 


where  I{r,z,r)  is  the  intensity,  =ho)  lcy„^p,  a„,p  is  a  cross  section  determined 


_  1  O  O 

empirically  [25,  42]  to  be  equal  to  6.4x10  cm  for  short  pulses  and  i  is  an  integer  denoting 
the  number  of  photons  needed  for  ionization,  i.e.,  i  =  lnt\lJiQy^  /  hco  +  1] .  The  characteristic 
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multi-photon  ionization  intensity  for  a  T  =  1/rm  laser  pulse  is  Ifj,p  =  5.8x10  W/cm  .  In  the 


tunneling  regime,  the  time-averaged  ionization  rate  for  a  linearly  polarized  laser  pulse  is  [40] 


^tun  ~ 


7/4,,  n1/4 


2  ^  ion 


,  (19) 


where  =  4.1x10^^  sec'*  is  the  fundamental  atomic  frequency,  /h  =  3.6x10*^  W/cm^  and  C/h 


13.6  eV  is  the  ionization  energy  of  hydrogen.  Finally,  since  the  peak  laser  intensity  in  a  guided 


1  "5  0 

filament  is  typically  ~10  W/cm  ,  the  ionization  process  is  neither  purely  multi-photon 
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(/<10^^  W/cm^)  nor  tunneling  (/>10^^  W/cm^).  In  the  intermediate  (/x  ~  1)  regime  an 

analytical  fit  is  employed,  having  the  form  Wx=cc^  I  ^«(-^))  jjj  expressions  for  the 

ionization  rates  ,  ajf,p  ,ax,T<\,  and  K2  are  fitting  constants  chosen  to  match  experimental 


measurements.  As  an  example.  Fig.  2  is  a  plot  of  the  ionization  rate  versus  intensity  for  a  laser 
wavelength  of  0.8  pm.  The  fitting  constants  for  this  plot  are  chosen  to  reproduce  the 
experimentally  measured  short-pulse  ionization  rate  reported  in  Refs.  25  and  42. 

ii)  Filamentation  in  Neutral  Air 

Perturbations  or  hot  spots  on  the  intensity  profile  of  a  laser  beam  can  grow  as  a  result  of  a 
filamentation  instability.  Filamentation,  i.e.,  transverse  break-up,  of  a  laser  beam  is  due  to  the 
interplay  between  diffraction  and  nonlinear  self-focusing.  Consider  a  laser  beam  propagating  in 
a  neutral  gas  for  which  the  nonlinear  focusing  power  is  Pnl  and  the  transverse  laser  intensity 
profile  is  slightly  perturbed  by  a  small,  localized  hot  spot.  The  spatial  growth  rate  of  this 
perturbation  due  to  the  filamentation  instability  [36,  43]  is  given  by 


r 

X  Inl  1 
tvxx^Pxl  x]_  j 


(20) 


where  I  is  the  laser  intensity  and  Xj^  is  the  characteristic  transverse  dimension  of  the  filament, 
i.e.,  spot  size.  As  a  function  of  the  dimension  of  the  filament,  the  growth  rate  vanishes  for 

1/9 

-  ^min  =  [FVz, ^  -^max  “  3^//(2FyVI ) 


xj^  =  V2  Xjjjin ,  and  decreases  inversely  with  x^^  as  Xj^  ->  qo  .  At  maximum  growth  rate  the 
power  within  the  filament  is  roughly  equal  to  Pnl-  It  is  therefore  expected  that  a  laser  beam  with 
a  power  P  will  break-up  into  filaments  where  N  <  P I  P^i  ■  As  an  example,  the  nonlinear 
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focusing  power  associated  with  air  for  a  1pm  wavelength  laser  is  in  the  range  of  -  2  GW; 
therefore  a  100  GW  laser  pulse  may  eventually  break-up  into  a  few  tens  of  filaments. 

iii)  White  Light  Generation 

The  nonlinear  interaction  of  an  intense,  short  laser  pulse  in  the  atmosphere  can  result  in 
significant  spectral  broadening  due  to  self  phase  modulation.  The  phase  of  the  laser  field 
becomes  modulated  through  the  time  dependent  refractive  index  by  nonlinear  effects,  ionization, 
Raman  processes,  etc.  Although  3D  effects  also  play  an  important  role  in  spectral  broadening  a 
ID  analysis  is  usefiil.  The  solution  of  the  laser  pulse  propagation  equation  in  ID  has  the  form 
A(z,t)  -  B(2,T)exp(i0(z,T))  where  the  amplitude.  Band  phase,  0  are  real  functions  of 
z  and  r .  The  instantaneous  frequency  of  the  pulse  on  axis  can  be  defined  as 


<y(z,  t)  =  <Uo 


d0(z,  t) 
dr 


(21) 


To  determine  dO!  dr  we  rewrite  the  full  nonlinear  propagation  equation  in  Eq.  (15)  in  the  form 


where  n{z,  r)  is  the  index  of  refraction  in  configuration  space  variables.  For  the  present 
purposes,  only  the  terms  contributing  to  the  index  of  refraction  from  bound  electrons,  Eq.  (7), 
Raman  scattering,  Eq.  (8),  and  the  plasma,  Eq.  (9),  will  be  retained.  The  bound  electron  and 
Raman  contributions  are  given  in  Eq.  (B15)  while  the  plasma  contribution  to  the  index  is 

-  o}^{z,r)l IHqCOq  .  The  nonlinear  index  can  be  written  as 


5n  —  w(z,  t)  YIq  —  Bl^hound  ^ Raman  ^plasma  ’  fhis  is. 


dn{Y,T)  =  n2l{r,T)  -  dT'W{r')R{r  -  r')/(r,r')  -  Q)^p{z,r)l In^col 


,  (23) 
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where  %  is  the  Raman  contribution  to  the  nonlinear  index  for  long  pulses,  W (r)  is  the 
population  inversion  variable. 


m  = 


^a>l+rf 


(Or 


e  sm{(ORr), 


(24) 


J 


00 


is  the  Green  function  for  the  Raman  process  normalized  such  that  jdTR(T)  =  I,  cor  is  the 

0 

characteristic  Raman  frequency  and  r2  is  a  phenomenological  damping  rate  (see  appendix  B  for 
the  details).  In  the  remainder  of  the  paper  we  will  assume  negligible  population  inversion, 

(r)  «  -1 .  To  obtain  the  instantaneous  frequency  along  the  pulse  we  substitute  the 
representation  of  the  complex  amplitude  A  into  Eq.  (22)  with  the  result 


oz  c 


OZ  c 


(25a) 


(25b) 


where  S (S rij )  is  the  real  (imaginary)  part  of  Sn,  which  is  assumed  small  compared  to  unity. 
The  instantaneous  frequency  spread  along  the  pulse  is  given  by  Eq.  (21),  together  with  Eq.  (25a), 


S  (o(z,  r)  =  (o(z,  t)-  cOq  - - 

r  J 


(Op  }  d5nr{z',T) 


dr 


(26) 


Note  that  only  the  nonlinear  terms  in  the  refractive  index  will  create  new  frequencies,  the  linear 
terms  redistribute  the  frequencies  within  the  pulse.  The  instantaneous  frequency  spread  is 
5(0  =  5(0f)pjf^({  +  5(ORpjfjpfj  +  5(0 ,  that  is 


5(t){z,  t)  = 


„5i?(r-r') 


1  d(o 


2  \ 


dr 


2n„a,l  Sr 


(OpZ 


(27) 


15 


where  I  =  cn^B^  l^rc  is  the  intensity. 

As  an  example,  consider  a  laser  pulse  with  wavelength  A  =  0.775  ^  w  ,  amplitude 
A(z,r)  =  -yJ^TT /g  / c  sm(n:T/Tj^),  for  0<r<r^,  and  zero  otherwise,  peak  intensity 

/^  =  5  X  lo'^  W/cm^  and  pulse  duration  =  500  fsec  propagating  in  air.  For  the  short  pulse 
regime,  the  bound  electron  and  Raman  effects  are  assumed  to  have  the  numerical  values 
n2=nR  =  3xl0"‘^cm^AV,  =  1.6x10^^  sec■^  and  r2  =  1.3x10^^  sec'*  [44],  For  the 

plasma  term,  the  dominant  ionization  mechanism  is  taken  to  be  multi-photon  ionization  of  O2  so 
that  the  ionization  rate  is  given  by  Eq.  (18)  with  f  =  8,  for  A  =  0.775 pm . 

For  these  parameters.  Fig.  3a  plots  the  individual  contributions  to  3i ,  given  by  Eq.  (23), 
due  to  bound  electrons,  Raman  scattering,  and  plasma.  Bound  electron  effects  produce  an 
increase  in  the  refractive  index  that  is  proportional  to  the  laser  intensity  while  the  generation  of 
plasma  causes  the  refractive  index  to  decrease  from  the  front  of  the  pulse  to  the  back.  The 
Raman  response  causes  an  increase  in  the  refractive  index  at  the  front  and  peak  of  the  pulse  and  a 
decrease  at  the  back.  The  sum  of  the  individual  contributions  to  the  refractive  index,  plotted  in 
Fig.  3b,  shows  that  for  these  parameters,  the  variation  of  the  refractive  index  is  of  the  order 

Figure  4a  plots  the  normalized  instantaneous  frequency  shifts  due  to  bound  electrons, 
stimulated  Raman  scattering,  and  plasma  after  propagating  for  50  cm  in  air.  For  these 
parameters  the  variation  in  the  frequency  shifts,  due  to  the  various  effects,  are  comparable  in 
magnitude.  The  bound  electrons  produce  a  red  shift  at  the  front  of  the  pulse  and  a  blue  shift  at 
the  back  while  ionization  produces  a  blue  shift  across  the  entire  pulse.  Stimulated  Raman 
scattering  produces  a  red  shift  near  the  front  of  the  pulse  and  a  blue  shift  at  the  back.  The  net 
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frequency  shift,  i.e.,  the  sum  of  the  bound  electron,  plasma,  and  Raman  contributions,  plotted  in 
Fig.  4b,  shows  a  20%  red  shift  at  the  front  of  the  pulse  and  a  larger  60%  blue  shift  at  the  back. 


IV.  SELF-GUIDED  PROPAGATION  OF  AN  IONIZING  LASER  PULSE 

In  this  section,  long  range  propagation  of  an  ionizing  laser  filament  in  air  is  considered. 
The  propagation  distance  of  a  laser  pulse  in  air  is  limited  by  a  number  of  processes.  Two 
fundamental  laser  pulse  propagation  mechanisms  that  can  in  principle  result  in  extended 
propagation  distances  are  i)  moving  foci  and  ii)  self-guiding.  In  the  moving  foci  mechanism  the 
focal  length  depends  on  laser  power  through  the  optical  Kerr  effect  and  different  temporal  slices 
of  the  laser  pulse  focus  at  different  distances  [45,46],  This  can  give  the  illusion  of  extended 
propagation.  However,  only  an  infinitesimal  fraction  of  the  laser  energy  is  propagated  over 
extended  distances.  In  the  self-guiding  mechanism  extended  propagation  distances  can  be 
obtained  by  balancing  the  defocusing  effects  of  diffraction  and  plasma  formation  against 
nonlinear  atomic  focusing,  i.e.,  Kerr  effect.  In  self-guiding,  losses  such  as  ionization  can  deplete 
the  laser  pulse  energy  and  significantly  limit  the  propagation  distance,  as  shown  below. 

The  process  of  ionization  and  optical  filament  propagation  can  be  analyzed  by  retaining 
diffraction,  the  nonlinear  refractive  index,  plasma  effects  and  laser  pulse  energy  depletion  due  to 
ionization  in  Eq.  (15),  i.e.. 


vi  +  2ik,^  +  y\A{r,z,rf  - 


Q)p{z,T) 


^Trih^  Ujpij  dj7g 


A(r,z,T)  =  0 


where  y 


n2  1 47r  c  . 


i)  Source  Dependent  Expansion  Method 
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The  following  analysis  of  Eq.  (28)  is  based  on  the  source-dependent  expansion  (SDE) 
method  which  was  originally  developed  in  Ref.  [47],  In  the  SDE  formulation  a  reduced  wave 
equation  of  the  general  form 


vi  +  2ik,^ 

oz 


A(r,z,T)  =  M(r,z,r)A(r,z,T), 


(29) 


is  solved  by  a  variation  of  parameter  technique  where  M(r,z,r)  is  a  known  nonlinear  function  of 


A(r,  z,  r) .  The  complex  electric  field  amplitude  is  given  by 


A(r,z,T)  =  B(z,t)  exp(i0(z,T))exp(-(l  +  ia(z,T))r^  /R^(z,t)),  (30) 


where  B  is  the  field  amplitude,  6>is  the  phase,  B  is  the  spot  size,  and  a  is  related  to  the  curvature 
of  the  wavefront.  The  quantities  B,  9,  R,  and  a  are  real  functions  of  z  and  r .  Using  the  SDE 
method  a  set  of  self-consistent  coupled  equations  for  the  pulse  amplitude  B(z,  r) ,  phase  9(z,  r) , 
curvature  a(z,  r) ,  and  spot  size  i?(z,  r) ,  can  be  derived.  Applying  the  SDE  method  we  find 


1  d(BR) 

BR  dz 

d6  (I  +  a^)  a  dR  1  da 
dz  R  dz  2  dz 


(31a) 

(31b) 


1  dR  2a 
Rdz  ^ 


-Gj, 


(31c) 


-  aGi,  (3 Id) 

2  dz  ko  R^ 

where  the  subscripts  (r,  i)  denote  the  real  and  imaginary  parts  of  the  function  respectively.  The 
details  of  a  related  derivation  using  the  SDE  method  can  be  found  in  Ref  [48].  The  complex 
functions  F  and  G  used  in  Eqs.  (3 1)  are  given  by 
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(32a) 


F{z,r)  =  ^\d{2r'^  IR^)M{r,z,T)txv{-2r^  IR^), 


G(z,t)  =  Y^jd(2r^  / R^)M(r,z,T)il  -  2r^  / ]exp(-2r^  /R^). 


(32b) 


Equations  (31)  can  be  combined  to  give  an  equation  for  the  pulse  power  and  spot  size, 

1  dP 


P  dz 


=  2R. 


(33a) 


^-“(l  +  KR^G,)+(2^  +  RGi 

dz^  klR^  '  \  dz 


Gi  +i?-^  =  0, 

dz 


(33b) 


where  P{z,r)-cR^  /16  is  the  laser  power.  For  the  present  problem  we  find  from  Eq.  (28) 


that 


M{r,z,r)  =  -  y\A{r,z,T)\  + 


2  a)p(r,z,T)  S/rik^  Ujon  dn^{r,z,T) 


Mf 


dz 


(34) 


Substituting  Eq.  (34)  into  Eqs.(32)  gives 


F(z,r)  = 


2k. 


-Lb^  +  _/  1  Ujondneq 


i  +  l 


c  i  B^  dT^ 


(35a) 


G{z,z)  = 


2k. 
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(35b) 


Finally,  substituting  Eqs.  (35)  into  Eqs.(33)  the  laser  pulse  power  and  spot  size  are  found  to  be 
given  by 
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(36b) 


where  Pj^i  -  c/2y  -  I  {271112)  is  the  nonlinear  focusing  power,  f  is  the  number  for  photons 

needed  for  multi-photon  ionization,  is  the  ionization  energy,  Wgo(z,r)  is  the  electron 
density  on  axis  generated  by  photo-ionization,  the  linear  index  has  been  set  equal  to  unity 
=  1 )  and  Tg  =  Imc^  -  2.8xl0~^^cm  is  the  classical  electron  radius.  The  on-axis 
electron  density  generated  by  multi-photon  ionization  is  given  by  Eq.  (25)  (evaluated  at  r  =  0) 


5t 


(f-l)! 


\^mp  J 


‘-n  5 


(37) 


where  I{z,r)  =  2P l{n]A?')  is  the  intensity  on  axis.  Equation  (36a)  indicates  that  the  laser 
power  decreases  as  a  function  of  propagation  distance  in  the  presence  of  ionization  since  energy 
is  expended  in  ionizing  the  air.  The  quantity  on  the  right  hand  side  of  the  equation  for  the  spot 
size,  Eq.  (36b),  represents  the  effects  of  energy  loss  due  to  ionization  and  indicates  that  a 
matched,  self-guided  filament,  is  not  possible.  The  single  photon  ionization  case  (f  =  1)  is  a 
special  situation.  The  right  hand  side  of  Eq.  (36b)  vanishes,  since  the  ionization  term  in  Eq.  (28) 
is  independent  of  r  for  f  =  1 ,  and  the  focusing  properties  of  the  pulse  are  not  affected  by  the 

ionization  process. 

ii)  Self-Guiding  Condition 

An  approximate  equilibrium  of  the  filament's  spot  size  can  be  found  when  the  energy 
depletion  due  to  ionization  is  low.  Neglecting  ionization  on  the  right  hand  side  of  Eq.  (36b)  and 
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taking  dR/ dz  =  0 ,  gives  an  approximate  equilibrium  condition  in  terms  of  the  filament  power, 
spot  size  and  on  axis  electron  density 

2?r  £  1 

Peq('r)  =  PnL  1  +  - "ectj)  ,  (38) 

I  V  +  lf  J 

where  the  subscript  eq  denotes  the  equilibrium  value.  The  electron  density  dependence  can  be 
removed  by  differentiating  Eq.  (38)  with  respect  to  r ,  and  using  Eq.  (37).  The  resulting 
equilibrium  condition  is  in  terms  of  a  differential  equation 


d  ( P(r)  -  P 

P(r)  j 


P\r) 
Rirf^  ’ 


(39) 


where  P(t)  =  Peq{'^)IPNL^  -^(0  =  Reqi'^) ! Peqi^)  and 


/3,„„  =  (0) 

((  +  1)2  (<-l)!  c2  /2 j 


where  =  {An  .  Since  there  is  no  plasma  at  the  head  of  the  pulse  to  counteract 


self-focusing,  the  initial  condition  on  the  power  is  P(0)  =  1 .  Consider  a  particular  set  of 
examples  for  which  the  spot  size  is  taken  to  vary  linearly  from  the  front  of  the  pulse  to  the  back, 
i.e.,  R{t)  =  1  +  er  Ixj^ .  Figure  5  plots  solutions  ofEq.  (39),  i.e.,  power  versus  r  for  several 
values  of  e  corresponding  to  cases  in  which  the  spot  size  is  increasing,  constant,  or  decreasing 
with  T .  The  laser  wavelength  is  taken  to  be  =  0.775  pm ,  Pnl  =  1  .7  GW ,  and  multi-photon 


ionization  of  O2  is  assumed,  i.e.,  f  =  8 .  To  counteract  plasma  defocusing  the  laser  power  must 
increase  with  r  to  maintain  equilibrium.  For  these  parameters,  the  variation  of  the  power  along 
the  pulse  is  the  smallest  for  the  constant  spot  size  example. 

iii)  Pulse  Energy  Depletion  Due  to  Ionization  and  Maximum  Propagation  Distance 
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The  rate  of  change  of  laser  pulse  energy  can  be  found  by  integrating  Eq.  (36a)  over  the 


pulse  length 


pulse  ^  nU ion 

az  2  f  J  ar  ’ 


(41) 


^i. 


where  E puise{z)  =  j6/rP(z,r)  is  the  pulse  energy  and  is  the  pulse  duration.  Equation  (41) 
0 

can  be  approximately  evaluated  by  taking  the  laser  spot  size  to  be  nearly  constant,  i.e., 
independent  of  r ,  and  using  the  approximate  equilibrium  condition  in  Eq.  (38), 


pulse  Uioyj  (£  +  1)" 


dz 


4r. 


PjTL) 

V  Enl 


-1 


(42) 


The  energy  loss  rate  is  independent  of  the  number  of  photons  needed  to  ionize  the  air  molecules 
when  i  is  large.  The  maximum  distance  a  pulse  can  propagate  Z-max  estimated  by 

assuming  that  all  the  pulse  energy  goes  into  ionizing  the  air. 


L 


4r„ 


E  pulse  (^) 


max 


Uio„(t  +  \)^{PIPNL-'iY 


(43) 


where  Epuise  (0)  is  the  initial  laser  pulse  energy.  In  principle,  extended  propagation  or  self- 

guiding  is  possible  only  when  the  laser  pulse  power  is  approximately  equal  to,  but  slightly 
greater  than,  the  nonlinear  focusing  power  which  for  air  is  2  GIF . 

iv)  Propagation  of  an  Ionizing  Laser  Pulse 

Equations  (36)  are  solved  numerically  to  illustrate  an  example  of  high  intensity  pulse 
propagation  in  air.  We  consider  a  laser  pulse  with  wavelength  X  =  0.775  pm  ,  and  a  uniform 

initial  spot  size  =  50  pm ,  which  corresponds  to  a  Rayleigh  length  of  =  ;r  /  A  «  1  cm  . 
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The  initial  intensity  at  the  front  of  the  pulse  is  (r  =  0)  =  4.3x1 0’^  W/cm^.  In  calculating  the 
ionization  rate,  we  assume  multi-photon  ionization  of  O2  so  that  f =8,  X-  0.775  pm  and 
ViQ„  =  12. 1  eV .  The  initial  power  profile  is  chosen  to  correspond  to  an  equilibrium  given  by  Eq. 
(38).  The  equilibrium  power  profile  is  perturbed  by  a  0. 1%  amplitude  modulation  and  its 
evolution  shown  in  Fig.  6.  For  the  first  10  Rayleigh  lengths,  power  is  depleted  uniformly 
throughout  the  pulse  by  ionization.  As  the  back  of  the  pulse  is  defocused  by  the  presence  of 
plasma,  the  intensity  drops  and  ionization  ceases.  The  front  of  the  pulse  remains  focused  and 
continues  to  ionize  and  lose  energy.  This  leads  to  a  localized  depletion  of  power  at  the  head  of 
the  pulse  which  is  evident  in  the  power  profile  at  z  =  40  Zr.  For  z  >  40  Zr,  the  pulse  intensity 
decreases  to  a  level  that  energy  losses  and  defocusing  due  to  ionization  become  negligible.  The 
result  is  that  the  power  profile  remains  relatively  unchanged  for  z  >  40  Zr.  The  spot  size  and 
intensity,  however,  continue  to  evolve. 

Figure  7  shows  the  evolution  of  the  laser  spot  size  at  three  different  locations  within  the 
pulse,  i.e.,  the  front  ( r  =  0 ),  within  the  body  of  the  pulse  (  r  =  1 2  fsec),  and  at  the  back 
(r  =  120  fsec).  The  spot  size  remains  relatively  constant  throughout  the  pulse  for  the  first  10 
Rayleigh  lengths  of  propagation  before  the  equilibrium  is  lost.  At  z  =  10  ,  the  front  and  back 

of  the  pulse  start  to  diffract.  The  spot  size  at  the  back  increases  at  a  faster  rate  due  to  the  initial 
defocusing  caused  by  the  plasma.  However,  the  portion  of  the  pulse  around  r  =  12  fsec 
maintains  at  a  constant  spot  size  for  over  200  Rayleigh  lengths  before  diffracting. 

Figure  8  shows  a  shaded  contour  plot  of  the  on-axis  laser  intensity  as  a  function  of  r  and 
z.  The  initial  laser  pulse  at  z  =  0  spans  the  length  of  the  plot.  The  defocusing  of  the  front  and 
back,  leads  the  formation  of  a  very  short  pulse  (~  5  fsec)  in  the  first  ~80  Rayleigh  lengths  of 
propagation.  However,  other  mechanisms,  which  are  not  included  in  the  reduced  SDE  equations. 
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e.g.,  group  velocity  dispersion,  will  affect  the  propagation  of  such  short  pulses.  This  short  pulse, 
which  is  characterized  by  P  ~  /’vl,  propagates  for  an  additional  150  Zr  before  diffracting.  This 
short  pulse  generation  process  resembles  the  relativistic  guiding  and  pulse  shortening  schemes 
proposed  in  earlier  works  [49,  50],  except  that  here,  ionization  contributes  to  the  defocusing  and 
guiding  is  accomplished  through  the  Kerr  effect  associated  with  bound  electrons. 

V  COMPRESSION  AND  FOCUSEVG  OF  LASER  PULSES  IN  AIR 

In  this  section  the  theoretical  model  together  with  simulations  are  used  to  study 
longitudinal  and  transverse  compression  and  ionization  of  chirped  laser  pulses  in  the  atmosphere. 
A  low-intensity  chirped  laser  pulse  propagating  in  air  can  compress  longitudinally  due  to  linear 
group  velocity  dispersion  and  focus  transversely  due  to  nonlinear  effects.  For  optimally  chosen 
parameters,  the  longitudinal  and  transverse  and  focal  distances  can  be  made  to  coincide  resulting 
in  a  rapid  intensity  increase  and  ionization  near  the  focal  region. 

The  propagation  of  the  high  intensity  laser  pulse  near  the  focal  region  is  markedly 
different  from  its  propagation  far  from  focus,  where  the  intensity  is  low.  In  the  following 
subsections  we  consider  separately  the  low-intensity  propagation  regime,  where  ionization  and 
Raman  processes  are  not  important,  and  the  propagation  near  focus,  where  the  laser  intensity 
becomes  sufficiently  high  that  ionization  occurs.  In  the  low-intensity  regime  a  coupled  set  of 
equations  for  the  laser  spot  size  and  pulse  length  are  derived.  Numerical  solutions  of  the  coupled 
equations  are  compared  with  the  full  numerical  3D  simulation  in  the  low  intensity  regime. 
Propagation  near  focus,  in  the  high  intensity  regime,  is  examined  using  the  full  3D  numerical 
simulations. 
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The  simulations  model  the  propagation  of -100  fsec  pulses  with  X  =  0.115 /urn ,  and 
intensities  as  large  as  -  10*^  W/cm^.  The  dominant  plasma  generation  mechanism  is  the  multi¬ 
photon  ionization  of  O2,  which  although  less  abundant  than  N2,  has  a  lower  ionization  energy. 
The  ionization  rate  given  by  Eq.  (18)  with  ^=8  agrees  well  with  the  ionization  rate  plotted  in 
Fig.  2  over  the  intensity  range  of  10*^  -  lO^'*  W/cm^.  The  parameters  used  in  modeling  the 
rotational  Raman  response  are  the  short  pulse  parameters  given  in  appendix  B,  i.e.,  rotational 

frequency  coji  =  16x10^^  sec'\  damping  rate  /’2  =  1.3x10^^  sec'^  and 

riR  K  ni  =  3xl0~^^cm^AV .  The  GVD  parameter  is  =2.2x10“^’  sec^/cm  and  higher  order 
dispersion  has  been  neglected.  For  the  range  of  intensities  and  pulse  durations  examined,  the 
effects  of  collisional  ionization,  recombination,  and  plasma  wakefields  are  not  important  and 
have  also  been  neglected.  For  a  typical  plasma  density  of  lO'^  cm'^,  peak  intensity  /-lO'"* 
W/cm^,  pulse  duration  -100  fsec,  and  laser  wavelength  X  =  0.775  pm,  the  ratio  of  the  pulse 
duration  to  the  plasma  period  is  tin  -  0.1  and  the  density  perturbation  associated  with  the 
Wakefield  is  of  the  order  dnjn^~\0'^'^.  A  brief  discussion  of  the  numerical  methods  used  in  the 

simulation  is  found  at  end  of  Sec.  II. 

i)  Low  Intensity  Propagation  Regime 

If  the  laser  intensity  is  sufficiently  low,  the  effects  of  ionization  and  Raman  scattering 
will  not  be  important  and  the  propagation  equation  in  a  spatially  varying  atmospheric  density, 
i.e.,  Eq.  (15),  reduces  to 

Vj_  +  2ik^-- - c^k„p2i^) — +  riz)\Air,z,T)f  A(r,z,r)  =  0.  (44) 
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In  writing  Eq.  (44)  higher  order  dispersion  has  been  neglected.  The  spatial  variation  of  the 
atmosphere  is  taken  into  account  mainly  through  the  z-dependence  of  the  group  velocity 
dispersion  coefficient,  ,  and  the  nonlinear  refractive  index,  n^{z),  both  of  which  are 
proportional  to  the  neutral  density.  In  principle,  the  wavenumber  is  also  dependent  on  the 

neutral  density,  but  is  taken  to  be  constant  since  the  fractional  change  in  is  less  than  -10“^, 
as  noted  in  the  discussion  of  Eq.  (3). 

Equation  (44)  can  be  solved  numerically;  however,  a  significant  simplification  of  the 
equation  is  possible  by  assuming  that  the  evolution  of  the  laser  pulse  is  self-similar.  That  is,  we 
assume  that  the  pulse  is  described  by  an  analytical  form  that  depends  on  certain  spatially 
dependent  parameters,  such  as  the  spot  size  and  pulse  duration  of  the  laser  pulse.  With  this 
assumption,  a  set  of  simplified  coupled  equations  can  be  derived  for  the  evolution  of  the  spot 
size,  pulse  duration,  amplitude  and  phase  of  the  laser  field.  Assuming  that  the  laser  pulse  has  a 
Gaussian  shape  in  both  the  transverse  and  longitudinal  directions,  the  complex  amplitude  can  be 
written  as 

=  B(z)  .  («) 

where  B  is  the  field  amplitude,  0  is  the  phase,  R  is  the  spot  size,  a  is  related  to  the  curvature  of 
the  wavefront,  T  is  the  laser  pulse  duration,  and  |3  is  the  chirp  parameter.  The  quantities  B,  6,  T, 
R,  a,  P  are  real  and  functions  of  the  propagation  distance  z.  The  instantaneous  frequency  spread 

along  the  pulse,  i.e.,  chirp,  is  Sco(z,t)  =  2/?(z)r/r^(z),  where  /?<0(>0)  results  in  a  negative 
(positive)  frequency  chirp,  i.e.,  frequency  decreases  (increases)  towards  the  back  of  the  pulse. 
The  full  frequency  chip  along  the  pulse  from  front  to  back,  i.e.,  from 
T  -  -T  to  T  -  T  ,\s  So)  full  -  A  pIT  . 
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Substituting  Eq.  (45)  into  Eq.  (44),  one  obtains  an  identity  in  the  variables  r  and  r . 
Equating  like  powers  of  r  and  t ,  the  following  coupled  equations  for  R  and  T  are  obtained. 


d'^R 


Eq  1 


52^ 

PNL^r 

d^T  _ 

1  ,  4/^2 

1  dP2 

dz^  ko  PnlR^T^  '  f" 


+ 


(46a) 


(46b) 


where  Eq  -  PT  is  proportional  to  the  laser  pulse  energy  and  is  independent  of  z. 


P{z)  -  nR{z)^  /(z)/2  is  the  laser  power,  I{z)  =  cnQB'^{z)l%K  is  the  intensity,  -PpL  ! ^ 


is  the  effective  self  focusing  power,  and  Pj^i  =1^  defines  the  usual  self-focusing 

power.  The  method  used  to  obtain  Eqs.  (46)  involved  equating  powers  of  rand  r ;  a  more 
rigorous  derivation  involving  the  source  dependent  expansion  method  [47]  would  result  in  P^^i 
being  equal  to  Pnl-  Hence,  in  numerical  solutions  of  Eqs.  (46)  we  set  =  Pj^i .  The  first 
term  on  the  right  hand  side  of  Eq,  (46a)  describes  vacuum  diffraction  while  the  second  term 
describes  nonlinear  self  focusing,  i.e.,  due  to  .  Nonlinear  self  focusing  dominates  diffraction 

when  P  >  .  The  curvature  parameter,  chirp,  phase  and  energy  evolve  according  to 


_  MM, 
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dz 
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respectively.  Equation  (47d)  indicates  that  the  pulse  energy  is  conserved. 

The  general  expression  for  the  Fourier  transform  of  the  pulse  amplitude  is 

00 

A(r,  z,a))  =  J dr  A{r,  z,  r) exp(ztyz-) .  Using  Eq.  (45)  we  find  that  the  Fourier  pulse  spectrum  is 


given  by 


S{co)  ~  exp 


4(1+/?^  (2)) 


The  full  Fourier  spectral  width  at  1/e  of  the  amplitude  is 


A£o(z) 


4(1  +  jp-{z)f'^ 


where  p{z')  is  given  by  Eq.  (47b).  Using  Eq.  (47b)  the  instantaneous  frequency  spread. 


8(o{z^r')  -  2^{z)r  IT^  {z) ,  becomes 


5o){z,  r) 


T  dT(z) 
dz 


Thus  the  frequency  chirp  vanishes  along  the  entire  pulse  when  the  pulse  length  reaches  a 


mimmum. 


For  a  spatially  uniform  neutral  density  and  in  the  absence  of  nonlinear  effects,  i.e., 
«2  =  0,  the  solution  ofEq.  (46b)  is  [38] 


V  /  n2 

r(z)  =  r^  ^ 

I  Zy-i  \^Lj ) 


and  the  chirp  parameter  is 


p{z)^ 


ip^dz 
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where  Tg  =  7(0),  Pq  =  ^{0),  and  Zj  =  1 2\(32  \  is  the  group  velocity  dispersion  length. 

Substituting  Eqs.  (51)  and  (52)  into  Eq.  (49)  yields  Aco  =  A{\  +  Po)  ITq,  which  indicates 
that  for  a  linear  atmosphere  the  Fourier  spectral  width  is  constant.  In  the  absence  of  nonlinear 
effects,  a  minimum  pulse  length  of  Tg  /(I  +  is  reached  eA.  z!  Z^.  =  -pQ  /(I  +  Po)’  provided 

the  chirp  parameter  is  initially  negative,  Po<^  ■ 

In  the  limit  that  the  pulse  length  does  not  change  appreciably,  it  can  be  shown  from  Eq. 
(46a)  that  the  variation  of  the  spot  size  with  propagation  distance  is  given  by 

(  p  V 

R(z)  =  R,i-2a„^+  ~  ,  (53) 

where  =  k^Rl  12.  Equation  (53)  shows  that  whenP  >  ,  the  spot  size  goes  to  zero  in  a 


distance 


(54) 


where  the  +  sign  is  chosen  so  that  z  is  positive. 

We  compare  the  solutions  of  Eqs.  (46)  and  (47)  with  the  full  numerical  simulation.  In  the 
present  example,  propagation  through  a  uniform  air  density  is  considered.  The  initial  laser  pulse 
at  z  =  0  is  described  by  Eq.  (45)  with  Z  =  0.775pm,  Rq  =lcm,  Tq  =0.66psec,  Pq  =-20  (a 


negative  chirp  with  \^5o) faiil cOg]^  »  0.05),  =  0  (a  collimated  pulse),  and  initial  peak  intensity 

lo  =  10^  W/cm^.  In  numerically  solving  Eqs.  (46),  the  nonlinear  index  of  air  is  taken  to  be 
-  6  X 10"’®  cm^AV,  which  is  larger  than  the  experimentally  measured  value  of  3  x  10"’® 
cm^AV  [44].  This  was  done  in  order  to  approximate  the  inclusion  of  Raman  effects  into  the 
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nonlinear  refractive  index.  Raman  effects  can  increase  the  nonlinear  refractive  index  by  an 
amount  equal  to  the  bound  electron  contribution,  e.g.,  see  Fig.  4.  The  ratio  of  the  peak  power  to 
the  nonlinear  focusing  power  is  P/Pnl  =  0.94  initially.  For  these  parameters,  the  focal  distance 
for  both  longitudinal  and  transverse  compression  is  expected  to  be  ~  0.5  km  based  on  Eqs.  (46). 

Figures  9a  and  9b  show  the  evolution  of  the  laser  spot  and  pulse  duration  as  the  focal 
point  is  approached.  The  pulse  length  decreases  almost  linearly  with  z  by  a  factor  of-  10  from  z 
=  0  to  z  -  0.49  km.  The  spot  size  decreases  at  a  slower  rate  from  z  =  0  to  z  =  0.4  km  and  then 
falls  rapidly  over  the  last  0. 1  km  from  focus.  The  corresponding  evolution  of  the  peak  intensity 
is  shown  in  Fig.  9c.  The  intensity  increases  relatively  slowly  over  most  of  the  propagation  path, 
gaining  a  factor  of  -  10  over  a  distance  of  0.4  km,  and  then  increases  rapidly  by  a  factor  of  >30 
in  a  distance  of-  0. 1  km.  The  reduction  in  the  spot  size  due  to  nonlinear  self-focusing  from  z  = 
0.4  km  to  z  =  0.49  km  is  mainly  responsible  for  the  enhanced  intensity  gain  observed  near  the 
focal  point.  The  results  of  the  full  simulation,  denoted  by  the  points  in  Figs.  9a  and  9b  and  the 
dashed  curve  in  Fig.  9c,  are  in  good  agreement  with  the  solution  of  Eqs.  (47).  Figure  10  shows 
surface  plots  of  the  distribution  of  laser  intensity  with  r  and  transverse  coordinate  x  at  z  =  0  and 
near  focus  at  z  =  0.49  km  obtained  from  the  full  simulation.  The  laser  pulse  focuses  both 
longitudinally  and  transversely  such  that  the  peak  laser  intensity  at  z  =  0.49  km  is  -  3x10^^ 

W/cm^,  which  is  a  factor  of  300  larger  than  the  initial  intensity;  note  the  change  of  scale  between 
Figs  10a  and  10b. 

ii)  High  Intensity  Propagation  Regime 

When  the  laser  pulse  is  sufficiently  intense,  ~10'^  W/cm^  for  -100  fsec  pulses,  ionization 
processes,  plasma  defocusing,  and  Raman  scattering  effects  become  important.  In  this 
subsection,  we  simulate  propagation  of  the  laser  pulse  from  the  previous  example  through  the 


30 


region  of  plasma  generation.  The  output  of  the  previous  simulation  is  extrapolated  for  ~2  m 
using  Eqs.  (46)  to  numerically  solve  for  the  laser  spot  size,  duration,  and  peak  intensity  close  to 
the  ionization  region.  The  extrapolation  was  performed  because  of  the  computational  difficulty 
in  using  the  full  nonlinear  simulation  to  accurately  model  the  large  variation  of  the  spot  size  near 
the  focal  point.  The  extrapolated  results  are  then  used  to  initialize  the  Gaussian  laser  pulse  in 
this  high  intensity  simulation  which  begins  at  =  492  m  .  The  propagation  distance  Az  is 

measured  relative  to  z^ ,  i.e  ,  Az  =  z-z^ .  The  pulse  at  Az  =  0  is  described  by  Eq.  (45)  with 
R  =  0.25cm ,  T  -  0.12psec ,  p  -1.2  and  a  =1 .  The  initial  peak  intensity  I  =  8.8x10*^ 

W/cm^  is  below  the  intensity  at  which  ionization  effects  become  important.  The  initial  ratio  of 
the  peak  power  to  the  nonlinear  focusing  power,  P/Pnl  “  2.5. 

Figure  1 1  shows  the  evolution  of  the  peak  laser  intensity  and  peak  plasma  density  near 
the  focal  point.  From  Az  =  0  to  Az  =  5  cm  self-focusing  causes  the  peak  intensity  to  increase  by 
a  factor  of  ~  8.  As  the  laser  intensity  reaches  the  ionization  threshold,  a  plasma  channel  is 
formed  which  is  highly  localized  near  the  laser  pulse  axis.  The  radius  of  the  plasma  channel  is 
~20pm .  Formation  of  a  plasma  channel  counteracts  the  nonlinear  focusing  effect.  From  Az  = 

5  to  Az  =  20  cm,  the  peak  intensity  is  limited  to  <7x10*^  W/cm^  while  the  peak  plasma  density 
is  on  average  ~  7.5x10^^  cm'^. 

Figure  12  shows  the  evolution  of  the  laser  intensity  profile  from  Az  =  6.3  cm  to  Az  =  16.8 
cm.  The  generation  of  a  highly  localized  plasma  channel  at  Az  =  4  cm  causes  the  trailing  edge  of 
the  laser  pulse  to  defocus  and  the  pulse  to  shorten  on-axis  as  shown  in  panel  (a).  Earlier  parts  of 
the  pulse  remain  focused  due  to  the  absence  of  plasma.  As  the  on-axis  pulse  length  decreases, 
the  plasma  density  also  decreases  thereby  allowing  the  trailing  parts  of  the  pulse  where  P  >  Pnl 
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to  refocus.  The  refocusing  of  the  trailing  edge  leads  to  the  double-peaked  intensity  profile  seen 
in  panel  (b).  Subsequently,  the  intensity  of  the  leading  peak  decreases  due  to  diffraction  until  the 
intensity  of  the  trailing  peak  becomes  the  global  maximum  [panel  (c)].  At  Az  =  16.8  cm,  the 
trailing  peak  is  reconstituted  in  such  a  way  that  the  laser  intensity  profile  appears  similar  to  that 
at  Az  =  7.6  cm.  In  earlier  works  [24]  it  was  proposed  that  these  recurrences  underlie  the 
experimentally  observed  long  distance  propagation  of  intense  pulses  in  air. 

Since  laser  energy  is  lost  to  ionization  whenever  refocusing  and  plasma  generation  occur, 
the  number  of  recurrences  will  be  limited.  Refocusing  will  not  be  possible  when  the  amount  of 
energy  lost  to  ionization  is  sufficient  to  cause  the  pulse  power  to  become  less  than  the  nonlinear 
self-focusing  power.  Figure  13  shows  the  pulse  energy  as  a  function  of  propagation  distance. 

The  sudden  decrease  in  energy  at  Az  =  5  cm  corresponds  to  the  location  at  which  a  plasma 
density  of cm'^  is  generated.  From  Az  =  5  cm  to  Az  =  25  cm,  ~7  %  of  the  pulse  energy  is 
lost  to  ionization.  Figure  14  shows  the  on-axis  profile  of  the  laser  intensity  and  power 
(normalized  to  Pnl)  at  Az  =  16.5  cm.  The  power  profile  at  Az  =  16.5  cm  remains  relatively 
unchanged  from  its  initial  Gaussian  profile  indicating  that  there  is  little  longitudinal  energy 
transfer.  Hence,  the  distortions  in  the  laser  intensity  profile  are  caused  mostly  by  transverse 
focusing.  Note  that  a  sufficient  amount  of  power  (P  ~  2  Pnl)  is  present  in  the  trailing  edge  of  the 
pulse  to  allow  for  self-focusing.  For  this  particular  example  however,  it  is  not  possible  to 
numerically  simulate  propagation  beyond  Az  =  25  cm  due  to  the  large  amount  of  spectral 
broadening  that  occurs,  pausing  the  theoretical  model  to  become  invalid,  e.g.,  at  Az  =  25  cm, 
\Sco\~  OAco^ . 

Figure  15  shows  that  the  on-axis  Fourier  spectrum  of  the  laser  pulse  broadens  with 
propagation  distance.  The  asymmetry  of  the  spectrum,  i.e.,  the  more  prominent  red  shift,  is 
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associated  with  the  gradient  of  the  laser  intensity  becoming  larger  at  the  front  of  the  pulse  than  at 
the  back  [38].  At  Az  =  17.7  cm  the  blue-shifted  part  of  the  spectrum  is  sufficiently  broad  that  it 
spans  the  visible  spectrum  indicating  that  white  light  is  generated  near  the  focal  point. 
Modulations  in  the  spectrum,  which  are  prominent  at  Az  =  15.2  cm,  are  caused  by  a  self¬ 
interference  effect  which  can  be  understood  as  follows.  Initially,  every  axial  position  within  the 
chirped  laser  pulse  has  a  different  frequency  as  shown  in  Fig.  16a.  As  the  pulse  propagates.  Fig. 
16b  shows  that  distortions  in  the  laser  envelope  cause  different  axial  positions  of  the  pulse  to 
have  the  same  frequency.  These  positions  represent  waves  with  the  same  frequency  but  with 
different  phase  that  can  interfere  constructively  or  destructively  depending  on  the  relative  phase 
difference.  This  interference  results  in  multi-peak  structures  in  the  spectrum  [51]. 
iii)  Vertical  Propagation,  Compression  and  Focusing 

We  use  Eqs.  (46)  to  examine  propagation  in  a  spatially  varying  atmosphere.  The  spatial 

variation  of  air  density  is  given  by  n^{z)  -  rt^(0)exp(-z/Z,^) ,  where  n^(0)  =  2.7  x  lO'^  cm^  is 
the  neutral  density  at  sea  level,  and  =  8  km  is  the  characteristic  scale  for  the  upward  variation 


of  the  air  density.  The  initial  laser  pulse  is  characterized  by  =  28  cm ,  =5  psec ,  lo  = 

3.2x10®  W/cm^,  A,  =  1.06  pm,  =  -46  (a  negative  chirp  with  do) y-y/zj/ qjq  «  0.02),  and 

=  10  (a  focusing  beam).  The  dashed  curve  in  Fig.  17a  shows  the  altitude  variation  of  the  air 

density.  The  solid  curve  shows  the  evolution  of  the  peak  laser  intensity  (normalized  to  its  value 
at  z  =  0)  with  altitude.  For  these  parameters  the  focal  length,  i.e.,  the  distance  at  which  the  pulse 
duration  and  spot  size  are  simultaneously  minimized,  is  ~  21  km.  The  intensity  increases  by  a 

factor  of  2x1 0'^  at  focus.  Most  of  the  intensity  gain  occurs  within  5  km  of  the  focal  region 
where  the  density  of  air  is  relatively  low. 
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Figure  17b  shows  the  variation  of  laser  spot  size  and  pulse  duration  with  z.  From  z  =  0  to 
z  ~  20  km,  the  spot  size  decreases  linearly  with  z,  indicating  that  the  transverse  focusing  is 
mostly  linear,  i.e.,  due  to  the  wavefront  curvature  of  the  pulse.  Nonlinear  self-focusing  becomes 
dominant  over  the  final  2  km  of  propagation  and  causes  the  spot  size  to  decrease  more  rapidly. 

In  the  absence  of  ionization  and  other  higher  order  nonlinearities,  the  spot  size  collapses  to  zero 
at  z  ~  22  km.  The  pulse  duration  decreases  continually  due  to  group  velocity  dispersion  and  is 
reduced  by  a  factor  of-  1/20  in  the  focal  region. 


VL  SUMMARY 

In  this  paper  we  have  investigated  a  number  of  key  physical  processes  associated  with 
short,  intense  laser  pulses  propagating  in  the  atmosphere.  Some  of  the  potential  applications 
stem  from  the  possibility  of  creating  an  atmospheric  ‘lamp’  at  a  remote  location  with  spectral 
characteristics  that  are  similar  to  a  white  light  source  [4,  8,  9].  The  applications  range  from 
remote  sensing  and  ultraviolet  fluorescence  spectroscopy  to  electromagnetic  countermeasures, 
hyperspectral  imaging,  differential  absorption  spectroscopy  and  induced  atmospheric  electrical 
discharges  (artificial  lightning)  [10-18,  52]. 

Nonlinear  equations  have  been  derived  that  include  dispersion,  nonlinear  self-focusing, 
stimulated  molecular  rotational  Raman  scattering,  multi-photon  and  tunneling  ionization,  pulse 
energy  depletion  due  to  ionization,  atmospheric  nonuniformity,  relativistic  focusing,  and  plasma 
Wakefield  generation.  The  nonlinear  equations  have  been  used  to  analyze  a  number  of 
phenomena,  such  as  the  compression  and  focusing  of  chirped  laser  pulses,  laser  filamentation, 
and  white  light  generation. 
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The  variation  of  the  nonlinear  refractive  index  and  instantaneous  frequency  shift  along  an 
intense  laser  pulse  in  air  has  been  calculated.  The  analytical  results  have  been  obtained  by 
assuming  a  self-similar  form  for  the  laser  envelope  and  deriving  coupled  envelope  equations  for 
the  amplitude,  phase,  curvature  and  spot  size.  As  an  example,  for  a  ~500  fsec  laser  pulse,  the 
frequency  shifts  associated  with  bound  electrons,  ionization,  and  Raman  scattering  are  found  to 
be  comparable.  The  combined  effects  of  bound  electron  nonlinearities,  ionization,  and  Raman 
scattering  produce  a  red  shift  at  the  front  of  the  pulse  and  a  larger  blue  shift  at  the  back. 

Using  the  source  dependent  expansion  method  (SDE)  a  coupled  set  of  equations  for  the 
spot  size,  laser  power  and  electron  density  is  derived.  A  necessary  condition  for  an  approximate 
equilibrium  of  a  single  optical  filament  has  been  derived  assuming  that  the  ionization  rate  is  low. 
The  equilibrium  involves  a  balancing  between  nonlinear  self-focusing,  diffraction,  and  plasma 
defocusing  and  is  shown  to  require  a  specific  distribution  of  power  along  the  filament.  For  the 
approximate  equilibrium,  the  laser  power  at  the  head  of  the  pulse  equals  the  nonlinear  focusing 
power  and  increases  towards  the  back  of  the  pulse.  The  increase  in  power  is  necessary  to 
counteract  the  effect  of  plasma  defocusing.  When  the  laser  intensity  is  sufficiently  high  to  cause 
considerable  energy  depletion  due  to  ionization,  an  equilibrium  solution  no  longer  exists  and 
self-guided  propagation  is  not  possible.  Numerical  solutions  of  the  SDE  equations  show  that  a 
~100  fsec  laser  pulse  can  propagate  in  a  self-guided  mode  for  ~10  Rayleigh  lengths,  after  which 
plasma  defocusing  erodes  the  front  and  back  of  the  pulse.  This  defocusing  leads  to  the  formation 
of  a  very  short  ~  5  fsec  pulse  which  can  propagate  for  ~  100  Rayleigh  lengths. 

A  method  for  generating  a  remote  spark  in  the  atmosphere  has  been  proposed  and 
investigated.  This  method  utilizes  the  dispersive  and  nonlinear  properties  of  air  to  cause  a  low 
intensity  chirped  laser  pulse  to  compress  both  longitudinally  and  transversely.  For  optimally 
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chosen  parameters,  the  transverse  and  longitudinal  focal  lengths  can  be  made  to  spatially 
coincide  resulting  in  a  rapid  intensity  increase,  ionization,  and  white  light  generation  in  a 
localized  region  far  from  the  source. 

Solutions  of  the  envelope  equations  are  found  to  agree  well  with  full  3D  simulations  in 
the  low  intensity  propagation  regime  far  from  the  focal  region.  Near  the  focal  spot  where  the 
intensity  is  large  enough  to  ionize  air,  plasma  filaments  can  be  generated  in  a  region  ~  1  m  in 
extent  and  significant  white  light  generation  occurs. 
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Appendix  A:  Derivation  of  Linear  Source  Terms 

This  appendix  outlines  the  derivation  of  the  linear  source  term  given  by  Eq.  (4),  Avhich 
describes  effects  associated  with  the  linear  polarization  of  bound  electrons  [38],  The  linear 
source  term  due  to  bound  electrons  can  be  written  as 


S^{r,z,t)  =  4nc^^d^^Vi{r,t)/dt^  =  e^  /2  +  c.c.,  (Al) 

where  is  the  linear  polarization  field,  the  phase  is  i^(z,  t)  -  kg  z  -  cOot ,  kg  and  cOg  are, 
respectively,  the  wave  number  and  fi'equency  of  the  carrier  field  and  is  a  unit  vector  in  the 
direction  of  polarization.  The  relationship  between  the  Fourier  transforms  of  the  linear 
polarization  field  and  the  laser  electric  field  E  is  given  by 

Pil**,®)  =  Xi(r,co)  E(r,o),  (A2) 

where  P^  and  E  are  the  Fourier  transforms  of  P^  and  E  ^  respectively,  and  is  the 

frequency  dependent  linear  scalar  susceptibility  which  may  also  be  a  function  of  r.  The 
convention  for  the  Fourier  transform  pairs  used  here  is 


(i-.ffl)  =  -7L  ]r,(T,t)e“  dt . 

y/'ZjZ  „oo 


P,(r,0  =  -^jP.(r,co)e-^“'c/© 
V27t 


The  relationship  between  P^  and  E  is  given  by 


(A3) 

(A4) 


(A5) 
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which  in  terms  of  Fourier  transforms  results  in  Eq.  (A2).  The  polarization  field  in  Eq.  (A5)  can 
be  written  as 

Pl  (•■>  0  =  -7^  f  Xi  (^)  E(r,  t  -  t)  ,  (A6) 

\2%  0 

where,  because  of  causality,  Xl  ("f)  =  0  for  r  <  0,  i.e.,  the  polarization  field  at  time  t  is  due  to 
the  electric  field  prior  to  time  t.  The  electric  and  polarization  fields  are  represented  in  the  form 


E(r,/)  =  4r,/>'>(">')e^/2  +  c.c.. 

(A7) 

+  c.a. 

(A8) 

where  Air,t)  and  C(r,t)  denote  the  complex  amplitudes  of  the  electric  and  polarization  fields, 
respectively. 

To  obtain  the  propagation  equation  describing  the  evolution  of  A(rj)  it  is  necessary  to 
express  C(r,  i)  in  terms  v4(r,/) .  Substituting  the  representation  for  the  electric  and  polarization 
field,  Eqs.  (A7)  and  (A8),  into  Eq.  (A6),  multiplying  both  sides  by  exp(/n)o0and  equating  like 
time  scales  we  find  that 


1 

C(r,()  =  -j^fzL(T)A(r,t-r)e“^‘’'^dT. 


(A9) 


Since  Xl(t^)  is  localized  near  rsO,  the  integral  in  Eq.  (A9)  can  be  approximately  evaluated  by 
expanding  the  electric  field  envelope,  A(r,/-r)  about  t  for  small  values  of  r. 


C(r,0  = 


j  XL  ix) 


1  -  r - h 


5  r'  53 


dt  2  6  Qp 


+ 


A{r,t)e'^°^  dr .  (AlO) 


Noting  that  the  derivatives  of  xi^o)  with  respect  to  0)q  are  given  by 
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V2it ;  ao„ 


(AH) 


where 


XeW  =  |xi{<<>J«""'' 


(A  12) 


substitution  of  Eq.(A12)  into  Eq.(AlO)  leads  to 


C(r  /)  =  y  d^hMd"A(r,t) 

„,0  «I  a®"  a" 


(A13) 


The  linear  source  term  amplitude  S^(r,t)  must  also  be  expressed  in  terms  of  the  electric 
field  envelope  A{r,t) .  Following  the  same  procedure  used  to  obtain  C(rj)  in  terms  of  A{r,t) , 
we  find  that 


I 

c'  J  4^  i  dx^ 


^2  LZw  m  ’ 


c"  £\  3(0 


(A14) 


where  we  have  assumed  that  5Xi,('r)/5T  vanish  at  r  =  0.  In  obtaining  Eq.  (A14)  the 


following  identity  was  used 


1 

J  Ph-^ 


(A15) 


The  linear  source  term  (r,  t)  can  be  rewritten  as 


O  /  A  f  ^o\  e  r  \  -e  0 

y(r,0=  —  ^ 

\  ^  J  £=0 


(A16) 


where  the  unitless  coeflRcient  «£  (r)is  given  by 
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a,(r) 


•(47rXi(coJ©^). 


(A17) 


In  terms  of  the  conventional  dispersion  parameters  Pi  [38],  defined  by 


Pi  =  ^{cD)ldco^  ,  where  P((b)  =  (© /c  )[1  +  ~  (a)/c)rt(G)),  the  coefficients 


are  given  by 


ai  - 


dco, 


(A18) 


0 


e.g.  «!  =  2 


^  \  ^2  =  1  -  c'^[p\  +  PoPi )  and  a3  =  -c'^(aj^p\p2  + 


1 - PqP\ 

®0  j 


For  the  cases  of  interest  it  is  sufficient  to  use  the  approximation  P\  ^  c 


-1 
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Appendix  B:  Derivation  of  Nonlinear  Source  Terms 

In  this  Appendix  the  nonlinear  source  amplitude  given  by  Eq.  (6)  is  derived.  It  is 
convenient  to  write  the  nonlinear  source  term  as 

~  ^bound  "*■  ^Raman  ^plasma  ^wake  ^rel  ^ion 


=  5'„(r.(y''<«'5,/2  +  e.c.. 


(Bl) 


where  the  various  contributions,  Abound  >  ^  Raman  ■>  ^  plasmas  ^wake  ^  ^rel’  ^  ion  ,  represent, 


respectively,  contributions  from  the  third  order  nonlinearity  associated  with  bound  electrons, 
stimulated  Raman  scattering,  plasma  generation,  wakefields,  relativistic  motion  of  free  electrons, 
and  laser  energy  depletion  due  to  ionization.  The  various  source  terms  are  derived  in  the 
following  subsections. 

a)  Nonlinear  Polarization  of  Bound  Electrons 

In  general,  the  nonlinear  source  term  due  to  bound  electrons,  Kerr  effect,  is 


^  hound  ^ ’0  ~  ^  ^bound^’^^^  ^  bound 


i  y/izj) 


ex/2 


+  C.C., 


(B2) 


where  the  phase  is  =kgZ  -  co^t .  The  nonlinear  polarization  field  is  given  by  [34,  35,  38] 

fbomd{r,t)  =  Xnl(^  ■  E>,  E  (r,(),  (B3) 


where  Xnl  ^^e  scalar  third-order  susceptibility  of  the  neutral  gas  and  the  brackets  (  denote 

a  time  average.  In  this  approximation  the  third  harmonic  component  of  the  nonlinear 
polarization  field  is  neglected  and  the  nonlinear  response  is  assumed  to  be  instantaneous.  The 
refractive  index  is  the  sum  of  the  linear  and  nonlinear  contributions  and  in  the  absence  of 
relativistic  effects  is  expressible  as 

n(r,  co)  =  (r,  a) +  n2  I ,  (B4) 
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where  /  is  the  time  averaged  intensity  and  is  the  nonlinear  refractive  index  associated  with 


bound  electrons.  Equation  (B3)  can  be  expressed  in  terms  of  n-i,  i.e.. 


^bound  (r,/)  =  =  c«2|^(r,0|^  E(r,0, 


(B5) 


where  is  the  linear  index,  «2  =  / ^o^]znl  >  (E  ■  E)^  =  /Stv  is 

the  intensity  and  |  ^2!  |  «  ~  ^  has  been  assumed.  Substituting  Eq.  (B5)  into  (B2)  and  using 

the  envelope  representation  for  E(r,  t) ,  i.e.,  Eq.  (2a),  the  amplitude  of  the  nonlinear  source  term 


IS 


s bound  ^  \A{r,t) 

An  c 

b)  Stimulated  Non-Resonant  Raman  Scattering 

Stimulated  Raman  scattering  of  laser  pulses  propagating  through  air  has  been  studied 
extensively  [28,  53-63].  For  altitudes  below  100  km,  the  dominant  Raman  process  for  long 
pulses  (  ~  nsec)  is  due  to  scattering  from  N2  molecules  involving  the  S(8)  rotational  transition 
from  the  J=8  to  J=6  rotational  states,  while  the  molecule  remains  in  the  vibrational  ground  state 
[55].  For  a  linearly  polarized  laser  with  wavelength  1.05  pm,  experiments  using  long  (~  nsec) 
pulses  indicate  that  the  Raman  gain  coefficient  is  ~2.5  cm/TW  [28].  The  observed  Raman  shift 

for  the  S(8)  transition  is  75  cm’’  (Am  ~  14x10^^  sec“^)  while  the  characteristic  relaxation  time 
for  excited  states  is  typically  0. 1  nsec  at  sea  level.  A  number  of  experimental  studies  have 
employed  shorter  (-100  fsec)  laser  pulses  to  investigate  Raman  scattering.  In  particular,  the  gain 
coefficient  and  damping  rate  have  been  measured  and  found  to  be  different  from  those 
appropriate  for  longer  pulses  [44]. 


A{r,t) 


(B6) 
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The  theoretical  model  used  here  to  incorporate  the  effects  of  stimulated  rotational  Raman 
scattering  into  the  general  nonlinear  propagation  equation  is  based  on  the  standard  density  matrix 
formalism  [64]  with  an  envelope  representation  for  the  laser  electric  field,  i.e., 

E(r,/)  =  +  c.c.,  where  A{r,t)  is  the  complex  amplitude,  ej_  denotes  a  unit 

vector  in  the  polarization  direction  and  c.c.  denotes  the  complex  conjugate.  The  envelope 
representation  allows  for  the  generation  of  a  multi-wave  Raman  spectrum,  i.e.,  harmonics  of  the 
Stokes  and  anti-Stokes  sidebands,  as  well  as  broadening  of  the  individual  lines.  The  model  also 
describes  the  Raman  response  in  the  transient  regime  and  accounts  for  the  natural  damping  or 
relaxation  of  excited  states  and  saturation  due  to  the  population  depletion  of  the  ground  state. 

In  our  model,  the  molecular  scatterer  is  assumed  to  have  two  nearby  rotational 
eigenstates,  1  (the  ground  state)  and  2,  with  corresponding  energy  levels  and  W2,  as  well  as 
an  excited  state,  e.g.,  an  electronic  or  translational  state,  with  energy  »  W2  -Wi.  Defining 
^nm  ~  ’  wherB  JO j  is  the  eigenff equency  associated  with  state y,  it  is  assumed  that 

Q3,,Q32  »  0)^  »  ,  where  is  the  fundamental  rotational  frequency.  That  is,  the 

Raman  process  is  non-resonant,  and  so  state  3  is  not  populated.  Because  of  the  assumed 
frequency  ordering,  we  can  effectively  take  w  032  s  Q .  Furthermore,  because,  eigenstates 

are  assumed  to  possess  definite  parity,  direct  transitions  between  states  1  and  2  are  forbidden. 

Raman  Stokes  scattering  consists  of  an  upward  transition  fi’om  state  (1)  to  a  virtual  level 
associated  with  state  (3)  followed  by  a  transition  from  the  virtual  level  down  to  state  (2).  In  the 
process,  a  photon  with  frequency  ojg  =  o)q  -  is  emitted.  Raman  anti-Stokes  scattering 
consists  of  an  upward  transition  from  state  (2)  to  a  virtual  state  followed  by  a  transition  fi'om  the 
virtual  state  down  to  state  (1)  thereby  emitting  a  photon  of  fi-equency  =  coq  +  ■^^21-  Since 
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the  population  of  state  (2)  is  much  smaller  than  that  of  state  (1)  in  thermal  equilibrium,  the  anti- 
Stokes  lines  are  generally  much  weaker  than  the  Stokes  lines  [34,  35], 

Stimulated  Raman  scattering  is  associated  with  a  nonlinear  polarization  field 

=  PRan,an{^’  eJ2  +  c.c. ,  which  glvcs  Hse  to  a  source  term  in  the  propagation 

equation  for  the  laser  envelope.  The  nonlinear  source  term  is  expressed  in  terms  of  the  nonlinear 
polarization  as 


(••>  ^  K  ^2  +  C.C.  .  (B7) 

Based  on  the  three-level  model  discussed  above,  it  can  be  shown  using  a  standard  density  matrix 
formalism  [64]  that  the  nonlinear  polarization  can  be  represented  as 


PRan.a„  =  Jfl  2(0  0  , 


(B8) 


where  Q(i)  is  the  unitless  oscillator  function  which  is  determined  by  the  system  of  equations 


dt^ 


+  +  +  2T, 


dt 


= 


n 


■W(t) 


(B9) 


6W  ^  nl  \A{r,t)\ 
dt  0)„ 


—  +  ^2  Q 


dt 


(BIO) 


where  is  the  Rabi  frequency  associated  with  the  peak  laser  amplitude  A^^ ,  and 

//  is  the  dipole  transition  moment  matrix  element  associated  with  transitions  to  state  3 .  The 
quantities  T^  and  r2  are  phenomenological  damping  rates  which  have  been  included 
heuristically.  The  quantity  W  is  the  difference  between  the  normalized  population  densities  of 
states  2  and  1,  and  JVq  =W(t^-oo) .  For  a  medium  in  which  all  molecules  are  initially  in  the 


ground  state,  Wg  =  - 1 . 
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Assuming  a  slowly  varying  amplitude  for  the  polarization,  i.e., 
d^PRaman^^^^  I  «  I  PRaman  L  the  source  term  for  stimulated  Raman  scattering  is  given  by 


(o: 


SRan,an  =  -  ^  Xt  Q{r,  0  A{r,  t) 


The  solution  of  Eq.  (B9)  can  be  rewritten  in  the  form 


27c  XL 


where 


Rit) 


+  rj '' 

COr 


e  sm{(ORt), 


is  the  Green  function  for  the  Raman  process. 


riR  = 


ISn'^XL  ^R 

cnl  n  0)1+  r}  ’ 


(Bll) 


(B12) 


(B13a) 


(B13b) 


is  the  Raman  contribution  to  the  nonlinear  index  for  long  duration  pulses  (t>/2  )  and  R(t) 


is  normalized  such  that,  ^dtR{t)  =  1 .  If  only  the  Raman  source  term  is  retained  on  the  right- 

0 

hand  side  of  Eq.  (15),  it  can  be  shown  through  a  stability  analysis  [34,  35],  that  an  initial  Stokes 
perturbation  d4(z  =  0)  on  a  CW  pump  laser  beam  can  grow  exponentially,  i.e., 

SA(z)  =  SA(0)exp(g/p  z),  with  a  maximum  gain  coefficient 


g  = 


8^^  XL  ^0  1 


C  nl  c  nrj  ’ 


(B14) 


where  /„  is  the  pump  intensity.  Stokes,  anti-Stokes  coupling  will  reduce  the  gain  coefficient. 
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The  nonlinear  polarization  field  in  Eq.  (B8)  has  a  contribution  which  is  third  order  in  the 
field  amplitude  and  thus  can  contribute  to  the  nonlinear  refractive  index,  «2  ■  Including  the 
bound  electron  and  molecular  Raman  response,  the  total  nonlinear  refractive  index  is 

nHL(tJ)  =  n2l(r,t)-nR\‘  dt'W(t’)R(t .  (B15) 

J  — 00 

For  a  constant  amplitude  laser  pulse  with  duration  vi ,  the  field  amplitude  can  be  written  as 
A  -  Aq  (6>(r)  -  ©{z-Ti)),  where  x  =  t  -  zlc .  We  assume  that  there  is  negligible  population 
inversion,  IT'(r)  «  -1 .  With  these  approximations,  Eqs.  (B13)  and  (B14)  indicate  that  the  total 
nonlinear  index  within  the  pulse  i.e.,  0  <  r  <  ,  is 

(l  -  {co%{cojix)  +  (r2  /®;i)sin((a^2:))),  (B16) 

where  can  be  written  in  terms  of  the  gain, 

eg.  (B17) 

CO^  +  r'2 

In  the  long  pulse  limit  (  r  »  1  /  r2 ),  rij^  represents  the  effective  nonlinear  index  due  to  Raman 
effects  and  the  total  nonlinear  refractive  index  is  =  «2  +f^R-  For  pulses  short  compared 
with  the  characteristic  Raman  times  ( r  «  1  /  ,  r  «  1  /  Ej ),  the  nonlinear  refractive  index  is 

due  to  purely  the  bound  electron  response,  i.e.,  -  «2  • 

Experiments  suggest  the  Raman  response  is  a  sensitive  function  of  the  pulse  duration. 

For  a  long  (~nsec)  laser  with  wavelength  \jum  (co^  =1.9x10^^  sec  ^),  the  rotational  Raman 
response  is  dominated  by  the  S(8)  rotational  transition  from  J  =  8  to  J  =  6,  which  is  characterized 

by  -1.4x10^^  sec'”^  r2 -lO'”  sec”',  and  g~2.5  cm/TW  for  air  at  STP  [28,56].  For  these 
parameters,  the  nonlinear  refractive  index  due  to  rotational  Raman  processes  is 
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nji  =  5,6  X  cm^AV  .  Assuming  that  bound  electron  and  Raman  effects  are  the  dominant 

contributions  to  the  nonlinear  refractive  index,  the  empirically  determined  value  of 

Unl  ~  ^2  the  long  pulse  regime  is  ~  5.6x10'*^  cm^AV,  giving  nj^l n^i  «  0.1.  More 

recent  experiments,  which  propagate  much  shorter  -100  fsec  laser  pulses  with  wavelength 
A  =  0.8  pm  through  air,  suggest  that  the  effective  parameters  for  the  short  pulse  regime  are 


r:  1.6x10'^  sec'‘,  «n2  ~3xl0  cm^AV,and  1.3x10'^  sec  ‘  [44]  giving  an 


effective  gain  coefficient  of  g  «  0.025  cm/TW  . 
c)  Plasma,  Wakefield,  and  Relativistic  Source  Terms 

The  source  term  in  the  wave  equation  due  to  the  motion  of  free  electrons  is  given  by 

^  free  ~  I  c^)  Si  I  j  where  the  plasma  current  density  J  satisfies  the  equation 


dJ  , 

dt  471 


2  r 


1  + 


bn. 


E(r,0. 


(B18) 


In  Eq,  (B 1 8),  (Op  (r,  /)  =  {A^rq^n^  (r,  t) /  rnf  is  the  plasma  frequency,  (r,  t)  is  the  electron 


density,  t)  is  the  plasma  density  perturbation  due  to  wakefields,  and  v^is  the  collision 


12  —1 

frequency  of  electrons  with  neutrals  and  ions,  for  air  Vg  «  3x10  sec  [19,  24].  Equation 
(B 1 8)  is  valid  even  if  the  plasma  density  is  increasing  because  of  ionization.  It  can  be  shown 
that  the  electric  field  associated  with  Sn^  (r,  t)  satisfies  the  equation  [65] 


+  c"^Vx  Vx  +  coi{r,t) 


dV 


= 


4W7 


{B19) 


The  Wakefield  density  perturbation  is  then  obtained  from  Sn^  =  V  •  I  An  q.  Writing  the  laser 


electric  field  and  plasma  source  terms  as,  E(r,/)  =  A{r,t)  /  2  +  c.c.  and 
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S frge{r,z,t)  =  S +  c.c. ,  respectively,  Eq.  (B18)yields 


'  ^  /  £•,-  C...  \t'  .  .  ^ 


free 


(Op{Y,t) 


1  + 


Sn^  5m 


rio 


m 


CO. 


A(r,t). 


(B20) 


0  J 


In  writing  Eq.  (B20),  the  plasma  current  density  J  was  obtained  to  order  Vg!  cOq  «  1  ■  The 
plasma  frequency  in  Eq.  (B20)  contains  contributions  from  ionization  and  relativistic  electron 
motion.  In  order  to  write  out  each  contribution  explicitly,  the  electron  mass  is  written  as 
m  +  5m,  where  m  is  the  electron  rest  mass  and  is  the  modification  due  to  relativistic  motion 
of  the  electron  in  the  laser  field.  To  second  order  in  the  field  amplitude  the  fractional  change  in 


the  electron’s  mass  is 


5  m/  m  =  — 
4 


r 


ymCCOo  J 


\A(r.lf. 


(B21) 


In  writing  Eq.  (B21)  it  has  been  assumed  that  |^zl / | «  1,  i.e.,  the  weakly  relativistic  limit. 
The  magnitude  of  qA!  me cOq  is  often  referred  to  as  the  laser  strength  parameter.  For  a  linearly 

polarized  laser  beam,  |  qA/mccoQ  |  =  8.6x10  A[pm]  /^^^[W/cm^],  where  X  is  the 
wavelength  in  microns  and  I  is  the  intensity  in  W/cm^. 

Using  Eq.  (B21),  Eq.  (B20)  can  be  written  explicitly  as  S free  =  S plasma  +  ^rel  +  ^wake  = 

where 


^plasma 


(oi{rj) 


; 


A{r,t) 


S 


rel 


q\A{r,t)\ 


4c" 


\2 


mccoQ  j 


A(rj), 


(B22a) 

(B22b) 
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A(r,0 


(B22c) 


_(0pir,i)3n^ 
^wake  o 


n. 


where  it  has  been  assumed  that  5m  Im  «  1 ,  v^l  (O^  «\,  and  \dn^\l  «  1 . 

d)  Laser  Energy  Depletion  Due  to  Ionization 

To  derive  the  source  term  describing  the  depletion  of  the  laser  energy  due  to  ionization 
we  note  that  the  rate  of  change  of  the  total  field  energy  (electric  and  magnetic)  due  to  only 
ionization  is 

dWf 

=  U ion  f (X,  y,  Z,  r)da ,  (B23) 

oz 

where  PFy  is  the  total  field  energy,  Ujon  is  the  ionization  energy,  and  dcx  is  the  differential 

cross  sectional  area.  The  total  field  energy  can  be  written  as 

Wj  =-\da\d{cT){^-'E)^l4n,  (B24) 

where  (  denotes  a  time  average,  (E  •  E)^  =\^\  ^2  and  -  cdr  is  the  differential  in  the  z 
direction.  Substituting  Eq.  (B24)  into  (B23)  we  find  that 

=  (B25) 

OCT 

Equation  (B25)  accounts  for  field  energy  lost  due  to  ionization  and  indicates  that  an  additional 
source  term  given  by 

should  be  present  to  properly  account  for  energy  depletion. 
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Chirped  Laser  Pulse  Air 


Plasma  Filaments  Optical  Filaments 

Figure  1;  Schematic  showing  filamentation  of  a  chirped  laser  pulse  in  air. 
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Figure  2:  Dependence  of  ionization  rate  on  laser  intensity.  Dashed  lines 
delineate  approximate  tunneling  (jk  <  0.5)  and  multi-photon  {jk  >  5) 
ionization  regimes. 
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Figure  4;  (a)  Instantaneous  frequency  shift,  Eq.  (27),  due  to  bound  electrons, 
Raman  effects,  and  plasma  versus  pulse  time  t  =  t  —  z/wg  after  propagating 
0.5  m  in  air.  A  laser  pulse  with  wavelength  A  =  0.775 /im,  peak  intensity 
=  5  X  lO^^W/cm^,  and  pulse  duration  tl  =  500  fs  is  assumed,  (b)  Total 
frequency  shift  (solid  curve),  i.e.,  SuJbound  +  Raman  +  ^^piasma  and  laser 
intensity  (dashed  curve)  versus  r. 
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Figure  5:  Dependence  of  normalized  laser  power  on  pulse  time,  t  —  t  —  z/vg, 
for  the  equilibrium  described  by  Eq.  (39).  The  spot  size  variation  is  assumed 
to  have  the  form  R{r)/Req  =  1  +  er/rn  with  tr  =  120  fsec  and  Reg  =  50 /xm. 
The  laser  pulse  has  wavelength  X  =  0.775  pm.  Multi-photon  ionization  of  O2 
(£  =  8  for  A  =  0.775  pm)  and  P^l  =  1  .7GW  are  assumed. 
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Figure  6:  Solutions  of  the  power  and  spot  size  equations,  Eqs.  (36),  showing 
power  profiles  vs.  pulse  time  r  =  t-  z/vg  sX  different  propagation  distances. 
The  power  profile  at  2:  =  0  is  an  equilibrium  described  by  Eq.  (38).  For  the 
initial  laser  pulse,  A  =  0.775  pm,  Io{t  =  0)  =  4.3xlO^^W/cm^,  To  —  120fsec, 
Ro{r)  =  50  pm,  i.e.,  an  initially  uniform  spot  size,  and  Rayleigh  length  Zr  = 
ttRI/X  =  1  cm.  Multi-photon  ionization  of  O2  (^  =  8  for  A  =  0.775  pm)  and 
=  1.7  GW  are  assumed. 
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Figure  7:  Normalized  spot  size  vs.  propagation  distance  at  three  different 
times  within  the  pulse  (r  =  0,  12,  120  fsec),  for  the  same  parameters  as  in 
Fig.  6.  The  portion  of  the  pulse  around  r  =  12  fsec  maintains  a  relatively 
constant  spot  size  for  over  200  Rayleigh  lengths. 
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Figure  8:  Level  plot  of  intensity  vs.  pulse  time  t  =  t  —  z/vg  and  normalized 
propagation  distance  z/Zr  for  the  laser  pulse  of  Figs.  6  and  7.  A  very  short 
])ulse  of  length  Ri  5  fsec  with  power  P  ^  P^l  is  generated  at  z  fs  80  Zr  and 
propagates  for  ^  150  Zr  before  diffracting. 


Figure  9:  Variation  of  (a)  laser  spot  size,  (b)  pulse  duration,  and  (c)  peak 
intensity,  from  z  =  0  to  z  —  0.49  km,  for  a  laser  pulse  with  initial  values 
Ro  =  I  cm.  To  =  0.66  psec,  Po  =  —20,  and  cvo  =  0.  The  initial  peak  intensity 
is  lo  =  lO^W/cm^  {Po/Pnl  =  0.94).  Solid  curves  denote  solutions  of  the 
envelope  equations,  i.e.,  Eqs.  (46).  Points  in  panels  (a)  and  (b)  and  the 
dashed  curve  in  panel  (c)  denote  full  scale  simulation  results. 


Figure  10:  Surface  plots  of  the  laser  intensity  at  z  =  0  and  z  =  0.49  km 
showing  compression  and  focusing  corresponding  to  the  simulation  of  Fig.  9. 
Peak  laser  intensity  increases  by  a  factor  of  300. 


Figure  11:  Full  scale  simulation  results  showing  the  variation  of  peak  laser 
intensity  and  peak  plasma  density  with  propagation  distance  in  the  vicinity  of 
the  focal  region.  The  laser  pulse  at  Az  =  0,  is  characterized  by  R  =  0.25  mm, 
T  =  120fsec,  P  =  7.2,  a  =  2,  and  peak  intensity  /  =  8.8  x  10^^  W/cm^. 
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Figure  13:  Variation  of  laser  energy  with  propagation  distance  in  the  region 
of  plasma  formation  for  the  simulation  of  Fig.  11. 


Figure  14;  Normalized,  on-axis  profile  of  laser  intensity  {I /I o)  and  power 
{PIPnl)  at  propagation  distance  Az  =  16.5  cm  for  the  simulation  of  Fig.  11. 
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Figure  15:  Fourier  spectrum  of  laser  intensity  as  a  function  of  wavelength 
at  propagation  distances  Az  =  0,  15.2  cm,  and  17.7  cm  for  the  simulation  of 


Figure  16:  Intensity  profile  (dashed  curves)  and  normalized  instantaneous 
frequency  spread  {Suo/uJo)  of  the  laser  pulse  at  (a)  A2:  =  0,  and  (b)  Az 
15.2  cm  corresponding  to  the  spectra  of  Fig.  15. 
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Figure  17;  (a)  Dependence  of  air  density  (dashed  curve)  and  peak  laser 
intensity  (solid  curve)  on  z  (altitude),  (b)  Normalized  laser  spot  size  (dashed 
curve)  and  pulse  duration  (solid  curve)  as  a  function  of  altitude.  Evolution 
of  the  laser  spot  size,  pulse  duration,  and  peak  intensity  for  the  given  density 
profile  are  calculated  using  Eqs.  (46).  Initial  conditions  are  given  by  Ro  = 
28  cm,  To  =  5  psec,  R  =  3.2  x  10®  W/cm^,  A  =  1.06 /rm,  /So  =  -46  (a  negative 
chirp  with  \Su!fuii\/(^o  ~  0.02),  and  cvo  =  10. 
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